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experiment  in  the  compressed  propellant  confined  by  steel.  The  gages  recorded 
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second  derivative  ofithe  Lagrange  pressure-specific  volume  compression 
curve  is  negative  if  this^'p^essure  region.  The  DDT  experiments  perform¬ 
ed  in  Lexan  tubes  with  propelTat^t  containing  50%  HMX  were  unsuccessful 
as  Lagrange  gage  experiments  becailse  the  flow  was  not  one  dimensional.. 
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experiments  with  propellant  containing  50%  HMX  confined  in  steel  tubes 
are  recommended  for  any  such  Lagrange  gage  experiments  that  may  be 
undertaken  in  the  future. 
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I  INTRODUCTION  AND  OBJECTIVES 


Research  on  the  deflagration-to-detonation  transition  (DDT)  is 
important  to  the  Air  Force  because  many  of  the  energetic  propellants 
required  for  present  and  future  long-range  delivery  systems  are  readily 
detonable,  can  be  used  as  good  explosives,  and  as  such  are  capable  of 
undergoing  DDT.  Rockets  motors  burning  HMX-based  propellant  have,  for 
example,  been  destroyed  by  DDT  on  several  occasions.  Possible  steps 
leading  to  detonation  in  the  rocket  motor  are  fracture  of  propellant 
ahead  of  the  flame  and  the  subsequent  formation  of  shock  waves  produced 
by  the  increased  burning  rate  of  fractured  propellant.  However,  in  this 
important  case,  the  mechanism  of  DDT  and  the  conditions  for  its  initiation 
are  not  adequately  understood. 

The  long-range  objective  of  the  present  research  program  is  to 
develop  a  computational  capability  for  assessing  the  DDT  hazard  in 
explosively  filled  propellants.  The  program  is  based  on  the  concept  that 
the  physical  and  chemical  processes  involved  in  DDT  must  be  identified, 
understood,  and  incorporated  into  a  realistic  model  of  DDT  before  the 
objective  can  be  achieved.  Theoretical  and  experimental  work  was  under¬ 
taken  to  model  and  establish  conditions  for  the  onset  and  occurrence  of 
DDT  in  HMX-based  propellants.  The  goal  of  the  theoretical  work  is  to 
formulate  a  realistic  model  for  the  flow  produced  by  an  unsteady  flame 
propagating  in  propellant  contained  in  a  closed  tube.  The  goal  of  the 
experimental  work  is  to  perform  Lagrange  stress  gage  experiments  to 
quantify  states  attained  ahead  of  the  flame  front  in  burning  HMX  propellant, 
to  establish  conditions  for  the  onset  of  DDT,  and  to  provide  data  for  the 
modeling  effort. 


II  THEORETICAL  WORK 


Our  theoretical  study  was  concerned  with  DDT  in  a  burning  propellant 
contained  in  a  closed  tube.  A  phenomenological  approach  was  taken  to: 

•  Develop  a  physical  understanding  of  the  thermo-hydrodynamics 
of  the  flame,  its  propagation,  and  the  flow  it  produces. 

•  Develop  a  model  of  the  flow  for  assessing  the  probability 
of  the  onset  of  DDT. 

To  develop  a  physical  understanding  of  the  thermo-hydrodynamics  of  the 
flame,  its  propagation,  and  the  flow  it  produces,  we  chose  a  model  for  the 
flame,  physically  interpreted  the  equations  governing  flame  propagation, 
and  derived  the  equations  governing  shock  formation  in  the  flow  in  front 
of  an  accelerating  flame.  To  develop  a  model  of  the  flow  for  assessing 
the  probability  of  the  onset  of  DDT,  we  constructed  a  solution  for  the 
unsteady  flow  produced  by  an  accelerating  flame  and  used  its  properties 
to  establish  conditions  necessary  for  the  onset  of  DDT. 

A.  Flames  Treated  as  Reactive  Discontinuities 

It  is  convenient  first  to  choose  the  simplest  model  and  treat  the 
flame  as  a  reactive  discontinuity.  Implicit  in  this  choice  are  the 
assumptions  that  the  burning  reaction  is  fast  and  that  the  hydrodynamic 
flow  produced  by  an  unsteady  flame  is  not  appreciably  influenced  by  the 
detailed  structure  of  the  flame.  In  this  case,  states  connected  by  the 
burning  process  are  governed  by  the  Rankine-Hugoniot  Jump  conditions 
expressing  the  conservation  of  mass,  momentum,  and  energy  across  the 
reactive  discontinuity. 

Let  v  denote  specific  volume;  F,  flame  velocity;  u.  particle  velocity; 
p,  pressure;  and  h,  specific  enthalpy;  and  let  the  subscripts  +  and  - 
denote  the  initial  and  final  states,  respectively  connected  by  the  flame. 
The  states  across  the  discontinuity  are  related  by  the  Rankine-Hugoniot 
jump  conditions^ 
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(1) 


v_(F  -  u+)  =  v+(F  -  u_) 

(u_  -  u+)2  =  (p_  -  P+)  (v+  -  v_)  (2) 

2(h_  “  h+)  =  (P_  "  P+)  (v_  +  v+)  (3) 

Equation  (1)  expresses  the  conservation  of  mass,  Eq.  (2)  expresses  the 
conservation  of  mass  and  momentum,  and  Eq.  (3)  expresses  the  conservation 
of  mass,  momentum,  and  energy.  The  locus  of  states  attainable  from  a 
given  initial  state  (p  ,  v, ,  h.)  lie  on  a  curve  in  the  (p,v)  plane  obtained 

T  T  *T" 

by  eliminating  h  between  Eq.  (3)  and  the  h  =  h(p,v)  equation  of  state  of 
the  burnt  products.  This  curve,  called  the  Hugoniot  curve  centered  on 
(p+,  v  ),  has  two  branches  when  the  discontinuity  is  reactive. 

When  the  reaction  is  exothermic,  the  branch  where  p  >  p+  and 
v  <.  v+  is  called  the  detonation  branch,  and  the  branch  where  p_  p+  and 
v  >  v  is  called  the  deflagration  branch.  With  the  present  model  of  the 

"•  T* 

flame,  burnt  states  are  represented  by  points  on  the  deflagration  branch 
of  the  Hugoniot  curve.  Moreover,  the  maximum  flame  speed  is  defined  by 
the  Chapman-Jouguet  (CJ)  point  where  F  =  u_  +  c_  (c  denotes  sound  speed), 
the  flow  is  sonic,  and  the  Rayleigh  line  through  the  initial  state  (p+,  v+) , 


is  tangent  to  the  Hugoniot  curve. 

Equations  for  deflagration  CJ  parameters  will  be  derived,  and  then 
the  differential  equation  governing  the  propagation  of  a  reactive  dis¬ 
continuity  will  be  formulated.  Equation  of  state  information  must  be 
specified  to  compute  CJ  parameters  and  states  on  the  deflagration  Hugoniot 
curve.  Let  the  subscript  o  denote  the  standard  state,  and  let  the  super¬ 
scripts  x  and  p  denote  propellant  and  propellant  products,  respectively. 
The  initial  enthalpy  of  the  propellant  and  the  h  =  h^(p,v)  equation  of 
state  of  the  products  are  written  as 
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(5) 


x 


AHX  +  hX 
o  + 


V 


and 

hp 


o 


+ 


kpv 
k  -  1 


(6) 


where  HP  =  AHP  -  k  p  v  /(k  -  1).  AH  denotes  the  heat  of  formation,  and 
k  denotes  the  poly tropic  index  of  propellant  products.  The  combination 
of  Eqs.  (3)  and  (6)  gives  the  equation  for  the  Hugoniot  curve  centered  on 
(P+v+)  as 


k  (yv_  -  v+)  +  P+(v+  +  vj  =  2  q  (7) 

where  y  =  (k  +  1) / (k  -  1)  and  q  =  [(AHX  -  AHP)  +  hx(p  ,v  )  +  kpv/k-l]. 

O  O  TT  OO 

At  this  stage  it  is  convenient  to  neglect  the  kp^v^/Ck  -  1)  term  in 

q  and  introduce  the  nondimens ional  variables  p_/p  =  p,  v  /v  =  V, 

2  2  -  +  -  + 
q/p  v  =  Q  and  U  =  (u  -  u  )  /p  v. .  Equations  (2)  and  (7)  can  then  be 

written  as 


(P  -  1)  (1  -  V)  -  u2 

(8) 

p(yv  -  1)  +  (V  +  1)  =  2Q 

(9) 

2 

Using  Eq.  (4)  and  the  definition  of  sound  speed  c  =  kpv,  the  CJ  condition 
F  -  u  =  c  can  be  written  as 


(10) 


The  CJ  parameters  are  obtained  by  solving  Eqs.  (8)  -  (10)  and  Eq.  (1). 
The  combination  of  Eqs.  (8)  ~  (10)  leads  to  the  following  equations 
relating  U  and  V  to  Q, 
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yU  +  (y  +  1)  -  2Q  =  0 


(1  -  V)  +  2B(1  -  V) 


■M 


B  =  0 


(11) 

(12) 


where  B  =  [2Q/(y  +1)  -  1],  The  roots  of  Eqs.  (11)  and  (12)  give  the 
particle  velocities  and  volumes  at  the  CJ  points  on  the  detonation  and 
deflagration  branches  of  the  Hugoniot  curve  centered  at  (p+,  v+) .  Solving 
Eq.  (11)  and  converting  back  to  dimensional  variables  gives  the  equation 


u_  =  u+  +  -k' y  [2{(k2  -  l)q  -  k(k  +  l)p+v+)}] 


1/2 


(13, 


where  the  positive  sign  gives  the  particle  velocity  at  the  CJ  point  on 
the  detonation  branch  of  the  Hugoniot,  and  the  negative  sign  gives  the 
particle  velocity  at  the  CJ  point  on  the  deflagration  branch. 

Whereas  liberation  of  chemical  energy  in  a  detonation  produces  mass 
motion  directed  toward  the  reaction  front,  the  liberation  of  chemical 
energy  in  a  deflagration  produces  mass  motion  directed  away  from  the 
reaction  front.  Expanding  the  square  root  term  in  the  equation  for  the 
roots  of  Eq.  (12)  gives  the  following  equations  for  the  CJ  volumes  on 
the  detonation  (14)  and  deflagration  (15)  branches  as 


v 


and 


v+  k  +  1 


,  2(k  -  1)Q 


k  +  1 


(14) 


(15) 


Rearranging  Eq.  (1)  to  give  the  equation 


F  -  u 


(u_  -  u+) 


+  (V  -  1) 


(16) 


2  2 

and  eliminating  (u_  -  u+)  by  using  Eq.  (13)  with  A  ■  2[(k  -  l)q  - 
k(k  +  l)p+v+]  gives  the  equation  for  the  maximum  deflagration  velocity 
as 
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F  -  u  - - f -  -  (li 

(1  +  A  /kp+v+) 

The  equation  for  the  CJ  deflagration  pressure  then  follows  from  Eq.  (9) 


P  -  1 


a2V+ 

(k  +  1)  (1  +  A2/kp+v+) 


When  A  /kp.v,  »1,  Eq,  (17)  gives  the  maximum  deflagration  velocity  as 

*r  *+• 

t 

P  V 
M-  + 

F  '  “+  *  k  -T*  (1 

and  the  corresponding  equation  for  the  change  in  pressure  across  the 
discontinuity  follows  from  Eq.  (19)  as 


P+  "  P_  k  +  1 


Equations  (13),  (15),  and  (18)  show  that  the  reaction  in  a  CJ 
deflagration  wave  produces  a  decrease  in  particle  velocity  and  pressure 
but  an  increase  in  volume.  Other  states  on  the  deflagration  Hugoniot 
curve  must  be  considered,  however,  because  the  CJ  deflagration  represents 
a  hypothetical  idealized  case.  These  states  lie  on  the  weak  deflagration 
bra  -.ch  where  pCj  <  p  _<  p+,  and  where,  according  to  Jouguet's  rules,  the 
flow  of  reaction  products  is  subsonic  with  respect  to  the  reaction  front. 
Jouguet's  rules  for  the  flow  in  front  of  a  weak  deflagration  give  an 
insight  into  the  problem  of  modeling  unsteady  deflagration  and  DDT.  The 
flow  ahead  of  the  wave  is  influenced  by  the  wave  itself  because  the  flow 
ahead  of  the  wave  is  subsonic  with  respect  to  the  wave  front.  Consequently, 
an  accelerating  flame  produces  an  unsteady  compression  wave  ahead  of  itself 
and  the  initial  pressure  for  the  burning  process  increases  as  the  flow 
develops.  The  burnt  states  are  represented  by  points  on  the  family  of 
deflagration  Hugoniot  curves  centered  on  an  adiabatic  compression  curve 
of  unburut  material.  Moreover,  the  continued  acceleration  of  the  flame 
will  result  in  shock  formation  and  initiation  of  detonation  in  material 
ahead  of  the  flame  front. 


'■  i*  J  :  1 


\  '  ‘  ■'  .  */ 

\ 

The  differential  equation  governing  the  propagation  of  a  reactive 
discontinuity  will  now  be  derived  to  obtain  a  more  quantitative  description 
of  unsteady  deflagration  waves.  Let  t  and  x  denote  time  and  Eulerian 
distance,  respectively,  and  p  -  p  =  Ap  and  u  -  u  ■  Au,  When  the 

*r  —  t  ^ 

deflagration  propagates  as  a  discontinuity,  the  time  variations  of  p+ 
and  p  satisfy  the  equations 


Dt 


.  * 


3t, 


3x+ 


(21) 


3E 

Dt  3t 


+  F 


i£. 

3x_ 


(22) 


Where  the  subscripts  +  and  -  denote  that  the  partial  derivatives  are 
evaluated  at  the  top  and  bottom  of  the  discontinuity  (the  initial  and  final 
states).  Substracting  Eqs.  (21)  and  (22)  arter  the  partial  time  derivatives 
have  been  eliminated  using  the  identity 


dt 


(23) 


leads  to  the  equation 


iE_+(F_u)|£  _(f_u)|£ 
Dt  dt,  dt  k  V  3x,  K  -  3x 


(24) 


A  similar  procedure  leads  to  the  corresponding  equation  for  u 


DAu  du  du  .  >  3u  v  3u 

dT  “  dT  -  Tt  +  <F  "  u+>  sT  '  <F  -  W 
+  -  + 


(25) 


The  equation  governing  the  discontinuity  is  obtained  by  combining 
(24)  and  (25)  with  the  equations  of  motion,  which  are  written  as 


du  3p 

dt  3x 


(26) 
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dp  2  dp  2  3u 

at  <=  at  c  psS  <27> 

with  p  =  1/v  used  to  denote  the  density.  Multiplying  Eq.  (25)  by  • 

M  -  p+(F  -  u+)  =  p_(F  -  u_)  and  adding  the  resulting  equation  to  Eq.  (24) 
leads  to  the  equation 

PAp  *  DAu  r  2  .2.  3u  ,  r  2  .2.  3u  /no. 

m  +  mE - p+[c+  -  <F  -  V  1  aS.  +  p-[c-  -  <F  -  “->  1  al  (28) 


Equation  (28)  gives  hydrodynamic  conditions  associated  with 
deflagration  waves.  Consider  first  a  deflagration  wave  propagating  at 
constant  velocity  when  D(Ap)/At  =  D(Au)/Dt  =  0  and  both  terms  on  the  right- 
hand  side  of  Eq.  (28)  are  zero.  Necessary  conditions  for  this  case  are 
either  3u/3x+  =  3u/3x  =  0,  or  3u/3x+  =  0  and  c  =  F  -  u  ,  because  unreacted 
material  satisfies  the  condition  c^_  >  F  -  u+.  Solutions  involving  a  pre- 
compression  shock  followed  by  a  deflagration  wave  can  be  constructed  to 
satisfy  both  of  these  conditions.  Rejection  of  the  CJ  condition  on 
physical  grounds,  however,  gives  3u/3x+  =  3u/3x  =  0  as  the  necessary 

hydrodynamic  condition  for  steady  deflagration  waves.  Of  particular 
interest  is  the  flow  when  the  particle  velocity  of  burnt  material  is 
zero.  In  this  case  the  strength  of  the  shock  and  deflagration  waves  are 
such  that  the  increase  in  particle  velocity  produced  by  the  shock  is 
exactly  offset  by  the  decrease  in  particle  velocity  produced  by  the 
deflagration. 


B.  Shock  Formation  in  Front  of  an  Accelerating  Flame 

The  factors  governing  shock  formation  in  front  of  an  accelerating 
flame  were  considered  because  formation  of  a  shock  is  necessary  for  DDT. 

A  shock  is  necessary  for  the  onset  of  DDT  because  detonation  is  a  shock 
supported  by  chemical  reaction.  At  this  stage  in  this  report  we  will  use 
the  subscript  o  to  denote  the  state  in  virgin  material,  rather  than  the 
standard  state.  With  the  notation  already  introduced  for  the  variables, 
we  let  h  and  a  denote  Lagrange  distance  and  strjss  and  write  the  equations 
of  motion  governing  our  one  dimensional  flow  as 


9 


(s);  '•  (S), 
(®)>  "••(«), 
(«)„--  (4 


It  is  convenient  in  considering  shock  formation  to  rewrite 
Eqs.  (29)  and  (30)  as 

-V°(*0„'(£)u  =  C“ 

v°  Wh = (£)„■  °« 

Equations  (32)  and  (33)  are  also  useful  for  the  Lagrange  analysis. 
The  quantities  C  and  are,  respectively,  the  slopes  of  the  curves 
of  constant  particle  velocity  and  stress  in  the  (t,h)  plane.  In 
other  words,  the  points  of  constant  particle  velocity  in  a  wave 
propagate  at  and  the  points  of  constant  stress  propagate  at  C^. 

A  shock  will  be  formed  in  the  compressive  part  of  a  wave  when  the 
wave  profile  steepens  as  the  wave  propagates.  The  conditions  for 
compressive  wave  steepening,  (3Cu/3u)^  >  0  and  (30^/30)^  >  0,  are 
thus  kinematic  conditions  for  shock  formation.  These  conditions  are 
related  to  dynamic  material  properties  by  the  following  equations 

/ 9C  \  (32u/?v2). 

/ _ u  ]  _  h 

\3u  /h  =  "  Vo  (9u/9v)h 
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I 


/M 

\3a  /h  =  vo  (30/8u)h  (35) 

obtained  by  differentiating  Eq.  (32)  partially  with  respect  to  u  and 

differentiating  Eq,  (33)  partially  with  respect  to  a.  Because 

(3u/3v)^  <  0  and  (3a/3u)^  >  0  during  compression,  the  dynamic  material 

properties  associated  with  shock  formation  follow  readily  from  Eq.  (34) 

and  Eq.  (35)  as  (32u/3v2),  >  0  and  (32a/3u2),  >  0. 

n  h 

Properties  of  the  Lagrange  stress-strain  curve  associated  with  shock 
formation  are  obtained  from  the  equation 


(36) 


obtained  by  combining  Eqs.  (32)  and  (33).  Differentiation  of  Eq.  (36) 
gives  the  equation 


2  02a/3v2)h 

Vo  (3c/3v)h' 


=  v  (3C  /3u),  +  C  (3C  /3a), 
o  u  h  u  a  h 


(37) 


The  condition  that  the  Lagrange  (°,v)  relationship  must  satisfy  for 

2  2 

shock  formation  to  occur  in  compression  follows  as  (3  a/3v  )^  >  0 
because  (3a/3v)^  <0. 

We  will  now  consider  shock  formation  in  a  simple  wave  adjacent  to 

a  constant  state.  In  a  simple  wave,  the  stress  and  specific  volume,  are 

2 

functions  only  of  particle  velocity  CQ  =  C^,  and  the  curves  of  constant 
particle  velocity,  of  constant  stress,  and  of  constant  volume  are  straight 
lines.  Equation  (36)  can  then  be  written  as 
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i1 .  i 


04 


and  differentiated  to  give  the  equations 


vq  (32a/3v2) 
2~  (3a/3v), 


.  C  /'  3C  \  /  3C  \ 


2  2 

These  equations  again  show  that  the  condition  (3  o/3v  ),  >  0  must  be 

1  h 

satisfied  for  shock  formation.  In  classical  work,  simple  waves  are 
treated  in  terms  of  characteristic  curves.  The  differential  equation 
for  the  characteristic  curves  is 


dt  — 


where  the  sound  speed  c  is  defined  by  the  identity  (c/v)  =  -  (3p/3v) 

s 

and  where  s  denotes  entropy  and  p  denotes  pressure.  The  characteristics 
associated  with  the  plus  sign  are  called  C+  characteristics  and  those 
associated  with  the  minus  sign  are  called  C_  characteristics.  At  this 
stage  we  identify  the  stress  with  the  pressure,  set  a  -  p  accordingly, 
and  use  the  e  =  e(s,v)  equation  of  state  to  rewrite  Eq.  (31)  as 


GO.- 


It  then  follows  from  Eq.  (36),  Eq.  (41),  and  the  identity 


/&\  o  (*2.\  +  (&\  (  ±l\ 

Wh  Ws  Wvl^'h 
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that 


(43) 


Equation  (40)  for  the  simple  wave  can  then  be  written  as 


dh 

dt 


=  +  C 
—  u 


(44) 


and  it  follows  that  the  characteristics  coincide  with  the  curves  of 
constant  particle  velocity  and  are  straight  lines. 

We  will  now  consider  shock  formation  in  a  simple  wave  produced  by 
an  accelerating  piston  or  an  accelerating  flame.  A  shock  is  formed  where 
the  forward  facing  characteristics  emanating  from  the  piston  or  the  flame 
intersect  to  form  an  envelope.  Integrating  Eq.  (44)  gives  the  particle 
velocity  field  for  the  flow  produced  by  a  piston  as 


h  =  Cu[t  -  tp(u)]  (45) 

where  the  subscript  P  denotes  the  piston,  and  gives  the  particle  velocity 
field  for  the  flow  produced  by  the  flame  as 


h  =  hf(u)  +  Cu[t  -  tf (u) ] 


(46) 


where  the  subscript  f  denotes  the  flame  discontinuity. 

The  condition  that  the  derivative  (3h/9u)t  vanishes  on  the  envelope 
formed  by  the  characteristics  gives  the  parametric  representation  of  this 
envelope  as 


C 

Ce  =  tp*'u^  +  (dup/dt)  (dCu/du) 


(47) 
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and 


e  (dUp/dt) (dCWdu)  (48) 

2  2 

where  t^u  du)  =  -  v  (3  p/3v  )  /2(3p/3v)  .  Equations  (47)  and  (48)  show 

OSS 

that  the  C+  .haracteristics  emanating  from  an  accelerating  piston 

dup/dt  >  0  form  an  envelope  in  front  of  the  piston  1m  =  0  when 

2  2  * 

(3  p/3v  )  >0.  Let  the  subscript  i  denote  the  initial  condition  at 

s 

the  origin.  When  (dUp/dt)^  >  0,  the  shock  is  formed  at  the  point  t  =  tc» 
h  =  cQtc  on  the  C+  characteristic  through  the  origin;  t  is  given  by 
Eq.  (47)  as 


Cc  =  (dup/dt) 1(dCu/du)i 


The  parametric  representation  of  the  envelope  associated  with  an 
accelerating  flame  is  obtained  from  Eq.  (46)  as 


(C  -  dhf/dt) 

te  tf^u^  +  (du_/dt) (d(^/du) 


dUU  (C  -  dh  /dt) 

he  =  hf(u)  +  Cu  (du_/dt) (dCu/du)  (51) 

We  consider  a  flame  with  an  initial  velocity  (dh^/dt)^  =  0  but  with  an 
initial  acceleration  (du_/dt)^  >  0.  In  this  case,  as  with  the  accelerating 
piston,  the  shock  is  formed  at  the  point  t  =  tc>  h  =  cQtc  the  C+ 
characteristic  through  the  origin;  the  time  t£  is  given  by 


_ _ o _ 

c  °  (du_/dt)i(dCu/du)i 


from  Eq.  (50). 


It  is  clear,  in  contrast  to  the  case  of  the  accelerating  piston, 

that  there  is  a  class  of  accelerating  flames  that  is  associated  with 

shock  formation  and  a  class  that  is  not.  These  classes  are  separated 

by  the  accelerating  flame  that  overtakes  the  first  characteristic  at  the 

time  tc>  A  shock  is  formed  by  accelerating  flames  that  would  overtake' 

the  first  characteristic  after  tc,  but  a  shock  is  not  formed  by  accelerating 

flames  that  overtake  the  first  characteristic  before  t  . 

c 

C.  Model  Solution  for  Assessing  the  Probability  of  the  Onset  of  DDT 

Our  model  solution  is  based  on  the  assumption  that  the  following 
chain  of  events  is  responsible  for  DDT  in  burning  propellant  contained 
in  a  closed  tube 

•  The  formation  of  an  accelerating  flame  by 
an  unsteady  burning  process. 

•  The  formation  of  a  shock  wave  in  unburnt  material 
by  the  compressive  action  of  the  accelerating 
flame  in  he  deflagration  to  shock  transition. 

•  The  formation  of  a  detonation  wave  in  material 
ahead  of  the  flame  by  the  shock  to  detonation 
transition  SDT. 

In  this  case,  DDT  will  not  occur  unless  the  accelerating  flame  produces 
a  shock  wave,  and  the  properties  of  this  shock  satisfy  the  conditions 
required  for  the  occurrence  of  SDT.  The  condition  for  shock  formation 
is  thus  not  a  criticality  condition  for  DDT,  and  the  shock  must  be 
characterized  to  assess  the  probability  of  DDT. 

Our  analytical  treatment  of  DDT  is  based  on  the  premise  that  the 

dynamics  of  the  flow  behind  and  in  front  of  the  accelerating  flame  must 

be  considered  to  characterize  the  shock  formed  in  front  of  the  flame  and 

to  assess  the  probability  of  DDT.  It  is  therefore  necessary  to  model  the 

accelerating  flame  and  the  unsteady  flows  on  both  sides  of  it.  The  model 

for  these  unsteady  processes  used  here  is  based  on  the  model  formulated  by 
3 

Jones  to  investigate  the  DDT  process  in  gases.  We  treat  the  flame  as  a 
reactive  discontinuity,  we  treat  the  flow  in  front  of  the  flame  as  a  simple 
compression  wave,  and  we  assume  that  the  density  gradient  is  zero  in  the 
flow  behind  the  flame.  The  one-dimensional  differential  equations  govern¬ 
ing  this  type  of  flow  in  a  closed  tube  can  be  solved  when  relationships  for 
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the  burning  rate  and  the  acceleration  of  the  flame  are  known.  This 
treatment  is  phenomenological  but  supersedes  previous  thermo-hydrodynamic 
treatments  by  making  the  states  attained  across  the  flame  compatible 
with  the  differential  equations  governing  the  conservation  of  mass  and 
momentum  in  the  flow  behind  the  flame,  and  by  making  the  states  compatible 
with  the  rear-boundary  particle  velocity  condition  in  the  closed  tube. 

It  is  convenient  in  constructing  solutions  for  such  flows  to  extend 
the  notation  intrrduced  on  page  8.  We  will  still  use  f  as  a  superscript 
or  subscript  to  denote  the  flame  discontinuity,  but  we  now  use  the  subscripts 
+  and  -  to  denote  quantities  in  front  of  the  flame  discontinuity  and  behind 
it.  States  across  the  flame  are  then  related  by  the  Rankine-Hugoniot 
jump  conditions  written  as 

P+  “  V>  “  Pf  <F  ~  uf)  (53) 

f  f  f  ,  f.  ,  f  fv 

P+  -  P_  =  p+  (F  -  U+)  <u+  -  u J  (54) 

2(e^  “  e+)  =  (P+  +  P*)  ~  v^)  (55) 

Again,  we  let  c  denote  the  sound  speed  and  the  subscript  o  denote  the  initial 
state  in  the  propellant.  We  also  assume  that  the  propellant  is  governed  by 
a  Murnaghan  equation  with  an  index  n  =  3.  The  equations  for  p  ,  P,,e  , 

T"  *T*  T 

and  the  slope  of  a  forward  facing  C+  characteristic  (u+  +  c+)  in  simple 
isentropic  compression  wave  can  then  be  written  as 


.  <p+  ■  !’o)v+  . 

%  + - 2 -  + 


Vo  -  V 


<v+  -  v0) 


(58) 


u,  +  c ,=  c  +  2u 
+  +  o 


(59) 


We  denote  the  extent  and  heat  of  reaction  by  X  and  q,  assume  that  the 
reaction  products  are  polytropic  with  an  index  k  =  3,  and  write  the 
equation  of  state  for  material  behind  the  flame  as 


a  =  e  -  Xq  + 

—  n  ^ 


,  s  Xp  vp 

(pv)  ro  0 
2  "  2 


-  (1  -  X) 


V+(p-) 


+  (1  -  X) 


(v (p  )  -  v  ) 
i  ~  o 


(60) 


where  v+(p  )  is  written  to  denote  that  the  volume  of  the  propellant  is 
related  to  the  pressure  by  the  same  equation  on  both  side  of  the  flame, 
and  vP  denotes  the  specific  volume  of  the  reaction  products  in  their 
standard  state.  Combining  Eq.  (55),  (58),  and  (60)  allows  us  to  write 
the  extent  of  reaction  in  the  flame  as 

Xf  =  AEf/Qf  (61) 


where 

Qf  ■  Uq  +  PQ(vP  -  vQ)  +  (c^Pq  -  3po)(v+(pf)  -  vQ)]  (61a) 

and 

AEf  =  [ (p_  +  p+) (v_  -  v+)  +  (pv)_  -  (pv)+ 

+(c^Pq  -  3pQ) (v+(p_)  -  VQ)]f  (61b) 

with  f  written  as  a  superscript  outside  the  bracket  for  notational  convenience. 
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With  x  and  t  used  to  denote  Eulerian  distance  and  time,  we  set 

f  3 

F  =  dx  /dt.  In  contrast  to  Jones  we  assume  that  the  burning  rate  is 
related  to  the  particle  velocity  and  the  pressure  at  the  flame  front  by 
the  equation 

F  **  u^  +  A(p^  -  pQ)  (62) 

where  A  is  the  burning  rate  parameter.  The  assumption  that  the  pressure 
exponent  in  Eq.  (62)  is  1  is  not  necessary,  but  is  made  here  for  convenience. 

We  are  now  in  a  position  to  consider  the  flow  produced  in  a  closed 
tube  by  a  flame  that  is  treated  as  a  reactive  discontinuity.  We  will 
first  derive  equations  for  the  initial  conditions  attained  in  the  tube 
and  then  present  solutions  for  flows  produced  by  flames  with  different 
types  of  acceleration. 

1.  Initial  Conditions  Produced  by  the  Propagation  of  a  Flame  in 
a  Closed  Tube 

The  flow  produced  by  the  propagation  of  a  deflagration  wave  from 
a  rigid  wall  is  of  particular  interest  because  our  experimental  study  of 
DDT  is  designed  to  model  this  rear-boundary  condition.  We  use  r  as  a 
subscript  and  as  a  superscript  to  denote  the  rear-boundary  and  we  use 
the  subscript  i  to  denote  initial  values  of  quantities  at  time  t  =  0. 

The  rear-boundary  condition  for  the  particle  velocity  at  the  wall 
is  expressed  by  the  equation 


u*  =  0  (63) 

Because  (u^)^  =  0  and  (p^)^  =  PQ,  it  follows  from  Eq.  (62)  that  F^  =  0, 

In  this  case  (u^)^  =  (u^)^,  and  it  follows  from  Eq.  (63)  that  (u^)^  =  0. 

Differentiating  Eq.  (63)  twice  with  respect  to  time  leads  to  the  conditions 
r  o  r  2 

du  /dt  °  d  u  /dt  =  0  and  it  follows  from  the  identities 


du*  9u_  9u_ 

IT  "  3t~  +  U-  hT 


(64) 
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and 


dV 

32u_ 

a2 

3  u 

2 

a2 

9  u 

du  l 

3u 

— 

r  — 

l  9ii 

+  (uf) 

-1- 

* 

dt2  ' 

3t2 

r 

+  2  -  3t3x 

r 

2  + 

3x 

r 

dt 

3x 

r 

(65) 


that  (3u  /3t)  =  (32u  /3t^)  =  0,  and  thus  that  (3u  /3t)^  =  (32'u_/3t2)^  =  0. 
We  now  use  the  following  identities  for  the  time  derivatives  of  uf  along  the 
flame  path. 


Du^  3u  3u 

Dt  3t^  +  ^  3x^ 


(66) 


2  f 
D  u 


Dt 


+  2F 


9  u_  2  3  U-  dF  ^ 

3t3x  +  ^  +  dt  3x 

f  9x  f  f 


(67) 


to  derive  an  equation  for  the  initial  value  of  the  particle  velocity 
gradient  behind  the  flame.  It  follows  from  equation  (66)  that  (Du_/Dt)i  =  0 
because  (3u  / 3t) ^  =  0  from  Eq.  (64)  and  F  =  0  from  Eq.  (62).  In  this  case, 
Eq.  (67)  gives  the  equation  for  (3u  /Sx)^  as 


because  (32u  /3 t2),  =  0  from  Eq.  (65).  It  follows  from  Eq.  (68)  that 

2  f  2  ~  1 

(Du  /Dt  )  and  (3u  /3x)^  have  the  same  sign  because  dF/dt^  >  0  is  the 
initial  condition  for  the  flame  to  propagate  along  the  tube. 

We  now  use  the  conservation  of  mass  expressed  by  Eq.  (53)  to 
calculate  the  initial  density  behind  the  flame  and  show  that  the  flame 
produces  a  compression  wave  in  virgin  material  as  it  propagates  from  the 
wall.  Setting  F^  =  (u^)^  =  (u^)  «  0  and  (Du^/Dt^  =  0  in  the  equation 

obtained  by  differentiating  Eq.  (53)  with  respect  to  time  leads  to  the 
following  equations  for  the  initial  density  behind  the  flame; 
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(p->i  -  Po 


(DF/Dt)i  -  (Du*/Dt)±' 


(DF/t), 


(69a) 


and 


<Oi  =  P0 


(DF/Dupi  -  1 
(DF/Du^)i 


(69b) 


f  f 

Because  (DF/Du.) .  >  0  it  follows  from  Eq*  (69b)  that  (p  ),  <  p  ,  and 
*  i  £  i  o 

it  then  follows  from  Eq.  (69a)  that  (Du+/Dt)^  >  0.  It  then  follows  from 
the  identity 


Dt 


du 


+ 


dt 


+  (F  -  up 


f-!!+ 

9x- 


(70) 


that  (du^/dt) ,  =  (Du^/Dt)  >  0.  Thus  the  particle  velocity  at  the  flame 

T  1  *T*  X 

front  starts  to  increase  and  the  deflagration  3tarts  to  produce  a 
compressive  wave  ahead  of  itself. 


The  equation  for  the  initial  value  of  the  particle  velocity 
gradient  behind  the  flame  is  obtained  by  combining  Eq.  (68)  with  the 
equation  obtained  by  differentiating  Eq.  (53)  twice  with  respect  to  time, 
and  making  use  of  the  continuity  equation  and  the  equations  expressing 
the  rear-boundary  conditions.  After  some  manipulation  the  equation  for 
(9u  /9x)^  can  be  written  as 


where  p.  is  defined  by  the  equation  p  =  p  p  .  It  follows  frotn  Eq.  (71) 

*r  *r  0  “r 

that  (9u  /9x)  <  0  in  a  closed  tube  because  the  term  in  braces  >  0  and 

j-  "  i 

(Du+/Dt)^  >  0.  Thus  the  flow  starts  to  develop  with  negative  particle 
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velocity  behind  the  flame.  In  the  case  when  F  satisfied  equation  (62)  and 
p+  and  p_._  satisfy  Eq.  (56)  and  (57),  the  initial  value  of  the  particle 
velocity  gradient  behind  the  flame  is  given  by  the  equation 


\  2  (Du>t)i 

Ap  c  +1 

0  0 

/i"  3  Co 

Ap  c 
o  0 

(72) 


Equations  for  the  pressure  are  readily  obtained  from  Eq.  (54) 

written  as 

pF  -  pF  =  M  (u*  -  uf)  (73) 

•  f  f 

with  M  =  p  (F  -  u  ).  Differentiating  Eq.  (73)  and  imposing  the  rear- 

T  T* 

boundary  conditions  leads  to  the  equation 


It  follows  from  the  identities 


f 

+  (F  -  V  9x7 

f 


(75) 


+  (F  -  uf) 


(76) 


that  (Dpf/Dt) j  -  (dp+/dt)i  =  (DpVDt)^  =  (dp^/dt)^  Ihe  asentropic 
condition  dp,  =  p  c  du ,  can  then  be  written  at  t  3  0  as 


Dt^  (3oCo 


(77) 
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and  combined  with  Eq.  (69)  to  obtain  the  equation 


(78) 


which  shows  that  the  pressure  starts  to  increase  behind  and  in  front  of 

the  flame  as  the  flame  starts  to  propagate  down  the  tube.  The  combination 

2 

of  the  energy  equation,  dp ,/dt  =  -  pc  (3u,/3x),  evaluated  at  t  *  0  with 

>  t*  *r  "T 

Eq.  (78)  then  gives  the  equation  for  the  initial  value  of  the  particle 
velocity  gradient  in  front  of  the  flame  as 


(79) 


f 

Because  Pq  >  (p_)^,  Eq.  (79)  shows  that  the  initial  value  of  the  particle 
velocity  gradient  in  front  of  the  flame  is  also  negative. 

Differentiating  Eq.  (73)  twice  with  respect  to  time  and  imposing 
the  initial  and  rear-boundary  conditions  leads  to  the  equation 


Equation  (80)  shows  that  the  pressure  difference  across  the  flame 
initially  increases  as  the  flow  develops. 

f 

2.  Flames  with  u+  Governed  by  a  Power  Law 

The  solution  with  the  particle  velocity  at  the  flame  front 
governed  by  a  power  law  of  the  form 


with  a  and  m  considered  as  parameters  was  presented  at  the  Seventh 
Symposium  (International)  on  Detonation  held  from  June  16-19,  1981 
at  the  United  States  Naval  Academy,  Annapolis,  Maryland.  Only  a  brief 
account  of  this  solution  is  given  here  (the  paper  presented  at  the 
symposium  is  included  as  the  Appendix  at  the  end  of  this  report) . 

Expressions  for  the  burning  rate  and  acceleration  of  the  flame 
given  by  Eqs.  (62)  and  (81)  were  used  to  construct  a  model  solution  for 
the  flow  produced  by  an  unsteady  flame  propagating  in  a  closed  tube. 
Equations  for  the  time  the  shock  is  formed  in  front  of  the  flame  were 
derived,  and  a  graphical  method^,  based  on  the  particle  velocity  pro¬ 
files,  was  used  to  determine  the  strength  and  position  of  the  shock 
at  later  times.  The  Hugoniot  equation,  Eq.  (55)  and  the  equation  of 
state  Eq.  (60)  were  used  to  calculate  the  time  t^  taken  for  the  reaction 

to  procede  to  completion  in  the  flame  discontinuity.  The  time  t  defines 

K 

the  time  interval  0  <  t  <  t  over  which  Eq.  (81)  can  be  used  to  describe 
the  acceleration  of  the  flame.  Consideration  of  the  evolution  of  the  shock 
preceding  the  flame  after  the  time  t  let  to  the  postulate  that  the  pressure 
attained  at  the  flame  front  when  the  reaction  proceeds  to  completion  in  the 
flame  can  be  used  to  assess  the  probability  <:f  the  occurrence  of  DDT. 
Examination  of  particular  solutions  in  the  light  of  this  postulate  led  to 
the  conclusion  that  DDT  depends  strongly  on  the  burning  rate  and  will 
probably  not  occur  unless  the  burning  rate  parameter  A  exceeds  a  critical 
value. 


3.  Flames  Attaining  a  Constant  Velocity  F 

For  the  sake  of  tractability,  the  model  solution  for  flames 
attaining  a  constant  velocity  was  based  on  the  following  equation 
governing  the  acceleration  of  the  flame. 


dt 

du 


a  ,^f  f. 
E  ■  ~2  <u+  *  "+> 


(82) 


where  a  denotes  a  characteristic  time  and  u7  is  the  particle  velocity 

A  T  A 

at  the  flame  front  at  the  time  t  =  t  when  F  =  F.  The  equation  for  the 


variation  of  the  particle  velocity  u£  with  respect  to  time  is  readily 


obtained  as 


t  =  t  - 


with  t  =  (a/'2)(u  /c  )  ,  by  integrating  Eq.  (82)  and  setting  t  =  t 

.  f  ~f  +  ° 

when  u ,  =  u , . 

+  + 

Our  model  solutions  are  also  based  on  the  assumption  that  the 

A  A 

flame  attains  its  constant  velocity  F  at  the  time  t  =  t  ,  when  the 

K. 

reaction  is  complete  in  the  flame  discontinuity,  and  thereafter  propagates 

*  f  Af  f 

at  F.  In  this  case,  the  condition,  u  =  u  when  X  =1,  can  be  used  to 

a^  a  a  T  i 

determine  the  values  of  u^,  F,  and  t  associated  with  a  prescribed  value 
of  the  burning  rate  parameter  A. 

We  will  now  derive  the  equations  governing  the  flame  trajectory 
and  the  states  attained  behind  the  flame.  It  is  convenient  to  denote  the 
equations  in  the  paper  in  the  Appendix  by  the  letter  A.  It  follows  from 
Eqs.  (A16),  (A17),  and  (82)  that  the  differential  equation  governing  the 
flame  trajectory  can  be  written  as 


B  -ft) *’6) 


where  D  =  3  +  c  /B  and  B  «  Ac  p  /3.  Integration  of  the  terms  obtained 
o  oo 

by  expanding  the  right  hand  side  of  equation  (84)  allows  us  to  write  the 
equation  for  the  flame  trajectory  as 


xf  =  I1  “  I2 


with 


and 
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Setting  u+  -  uu  in  Eqs.  (85),  (86),  and  (87)  gives  the  distance  taken 
for  the  flame  to  attain  the  constant  velocity  F  as 


Again  the  calculation  of  states  behind  the  flame  is  based  on  the 
assumption  that  (3p_/3x)  =  0.  In  this  case,  the  mass  behind  the  flame 
can  be  written  simply  as  p_xf,  and  formulation  of  the  conservation  of 

mass  for  the  material  that  has  passed  through  the  flame  front  leads  to 

_  £ 

the  following  equation  for”p  x 


Because  p  =  p  ,  the  density  calculated  from  Eq.  (89)  can  be  used  to 

f  ”  f 

calculate  u_  and  p_  from  the  jump  conditions,  Eqs.  (53)  and  (54),  and 
these  quantities  can  be  used  to  complete  the  characterization  of  the 
u_  (x,t)  and  p_(x,t)  fields  behind  the  flame.  We  set  M(uf)  =  p^(F  -  uf) 

for  convenience,  and  use  Eqs.  (62),  (56),  and  (57)  to  obtain  the  following 

•  f  f 

expression  for  M(u.)  in  terms  of  u., 

T  *T 


M(up  = 


/f\T  /f\  /  fv2  /  f  ' 

■•■•ft)- ft) -ft)  *(? 


(90) 


It  is  convenient  after  performing  the  integration  in  Eq.  (89)  with 

f 

Eq.  (82)  and  (90)  to  write  the  equation  for  p  x  as 


p  x  =  G1  -  G2 


(91) 


with 


gi  =  pb«(~)(~)  l  +  2f“)  +  f~±>)  +  ^ 

1  °  \Co/\Co/  2  \Co/  VCo/  5lCo/ 
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(92) 


and 
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-Hk) 


(93) 


The  equation  for  the  density  p  when  the  flame  attains  the  constant 
velocity  F  is  readily  obtained  from  Eqs.  (91)-(93)  as 


Ba 

p-  xf  =  po  1 


/\  r  -*  f  ^  r  l-  /N  f 

(5) 


(94) 


and  combining  Eq.  (88)  and  Eq.  (94)  gives  the  corresponding  equation 
for  p  as 


P  =  P, 


1  +  (uf/c  )  +  2/5  (u^/c  )2  +  1/15  (uf/c  )3 
D/3  +  1/2  (u*/c  )  +  1/10  (ui/c  )2 

TO  TO 


(95) 
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We  are  now  in  a  position  to  construct  particular  model  solutions 
for  flames  attaining  a  constant  propagation  velocity.  Values  of  the  input 
parameters  used  to  generate  these  solutions  are  p  =  1.6  g/cc, 

n  o®  o 

c  =  2  mm/ys,  q  =  54.69  kbar  cc/g,  vp  =  8.4  x  10  cc/g,  a  =  10J  ys, 

p  =  1  x  10  kbar,  and  values  of  the  burning  rate  parameter  A  in  the 
0  -3  -2  . 

range  1  x  10  -  1.2  x  10  mm/ys  kbar.  Equations  (95),  A29),  (61),  (61a) 

f  Af  f 

and  (61b)  and  the  condition  u  =  u  ,  when  X  =1,  were  used  to  calculate 

/>.£  T  1" 

values  of  u  for  the  prescribed  values  of  the  burning  rate  parameter  A,  and 

Af 

the  corresponding  values  of  p  were  used  to  assess  the  probability  of  the 

+  f  f  Af 

occurrence  of  SDT.  For  values  of  u  in  the  range  0  <  u <  u, ,  the  flame 
trajectory  was  calculated  from  Eqs.  (85)-(87),  and  the  states  connected 
by  the  flame  were  calculated  from  Eqs.  (85) — ( 87 ) ,  (91)-(93),  (A28),(A29), 
(61a),  and  (61b).  Values  of  t,  u+,  u  ,  p+,  and  p  calculated  with 
different  values  of  the  burning  parameter  are  given  in  Table  1. 


Table  1 

A  /\  /-  A  r*  A  r.  A 

.  f  ±  r  f 

t,  U  ,  U_  ,  p  ,  p_ 

CALCULATED  FOR  VALUES  OF  THE  BURNING  RATE  PARAMETER  A 

A  Af  Af  Af  Af 


t 

u+ 

X. 

u_ 

P+ 

j- 

P_ 

mm/ys  kbar 

ys 

mm/ys 

mm/ys 

kbar 

kbar 

5.5  x  10~3 

10,1 

0.284 

-0.043 

10.4 

10.1 

1.0  x  10~2 

44.4 

0.596 

-0.205 

25.3 

21.1 

1.2  x  10“2 

80.2 

0.801 

-0.381 

44.4 

25.4 

Af 

Examination  of  Table  1  shows  that  the  values  of  the  shock  pressures  p, 

T 

for  assessing  the  probability  of  DDT  for  the  values  of  the  burning  rate 
-3  -2  -2 

parameter,  5.0  x  10  ,  1,0  x  10  ,  and  1.2  x  10  mm/ys  kbar,  are 

respectively  10,4,  25.3,  and  44.4  kbar.  DDT  is  expected  to  occur  when 

Af 

the  shock  pressure  p  produced  by  the  flame  is  44.4  kbar  and  perhaps 

Af 

when  p,  =  25.3  kbar,  but  not  when  p  =  10.4  kbar.  Here  again,  as  in  the 

+  +  Af 

calculation  presented  in  the  Appendix,  the  dependence  of  p+  on  A  shows 
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that  the  probability  of  DDT  depends  strongly  on  the  burning  rate  parameter, 
and  also  leads  to  the  conclusion  that  DDT  will  probably  not  occur  in  a 

_3 

closed  tube  when  A  <  0.9  x  10  mm/ys. 

Theoretical  Results  and  Conclusions 

Equations  governing  flames  treated  as  reactive  discontinuities  were 
derived  and  physically  interpreted. 

Conditions  for  shock  formation  were  established,  and  the  compressive 
flows  produced  by  an  accelerating  flame  and  by  an  accelerating  piston 
were  compared  to  provide  a  better  physical  understanding  of  the  compressive 
action  of  an  accelerating  flame. 

Solutions  describing  the  one-dimensional  unsteady  flow  produced  when 
an  accelerating  flame  in  a  propellant  or  explosive  contained  in  a  closed 
tube  attains  a  constant  propagation  velocity  were  constructed  to  model  DST 
and  gain  an  insight  into  DDT.  The  flame  was  treated  as  a  reactive 
discontinuity,  the  flow  in  front  of  the  flame  was  treated  as  a  simple 
compression  wave,  and  the  density  gradient  was  assumed  to  be  zero  in  the 
flow  behind  the  flame.  The  description  of  the  flow  is  phenomenological 
because  the  burning  rate,  the  acceleration  of  the  flame,  and  the  density 
gradient  behind  the  flame  must  be  prescribed  to  construct  the  solution. 

But  this  treatment  supersedes  previous  thermodynamic  treatments  by  account¬ 
ing  for  the  flow  behind  the  flame  as  well  as  for  the  flow  in  front  of  it. 

Particular  solutions  were  constructed  for  different  values  of  the  burning 

~f 

rate  parameter  A.  The  value  of  the  pressure  at  the  flame  front  p ,  when 

1  A 

the  flame  attained  its  constant  propagation  velocity  after  the  time  t 

was  used  to  assess  the  probability  of  the  occurrence  of  DDT.  The  dependence 
~f 

of  p+  on  A  shows  that  DDT  depends  strongly  on  the  burning  rate  and  will 

probably  not  occur  unless  the  burning  rate  parameter  exceeds  a  value  of 

_2 

approximately  1.0  x  10  mm/ys  kbar. 

The  model  solutions  presented  in  this  report  for  the  flow  associated 
with  an  accelerating  reactive  flame  discontinuity  provide  guidelines  for 
assessing  the  probability  of  the  occurrence  of  DDT  in  a  burning  propellant 
contained  in  a  closed  tube.  Questions  regarding  the  applicability  of 
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these  solutions  to  real  systems  must  be  considered,  however,  because  the 
extent  of  reaction  in  the  flame  discontinuity  increases  as  the  flow 
develops.  In  this  case,  the  global  assumptions  made  to  model  the  flow 
may  be  reasonable  for  composite  propellant  containing  components  with 
markedly  different  reaction  rates. 

A  more  fundamental  theoretical  study  of  the  flow  produced  by  an 
unsteady  flame  must  be  undertake  to  obtain  a  more  basic  understanding 
of  DDT.  Questions  regarding  the  treatment  of  a  flame  as  a  reactive 
discontinuity  and  the  application  of  a  steady-state  burning  rate 
expression  to  an  unsteady  flame  must  be  addressed  in  this  study,  and  a 
more  realistic  model  must  be  developed  for  predicting  the  trajectory  of 
the  unsteady  flame  and  the  flow  fields  it  produces. 


Ill  EXPERIMENTAL  WORK 


Initially  the  experimental  studies  were  delayed  because  we  could  not 
find  a  source  of  propellant  for  the  Lagrange  gage  experiments.  Protracted 
negotiations  to  obtain  HMX-based  propellant  for  the  DDT  experiments  from 
Lawrence  Livermore  National  Laboratory  LLNL  were  unsuccessful  for  the 
following  reasons: 

•  Propellant  was  in  short  supply. 

•  The  machining  costs  and  times  to  meet  tolerances 
in  the  dimensions  of  propellant  samples  (needed  to 
construct  the  DDT  gage  assemblies)  were  prohibitive. 

We  subsequently  found  a  source  of  HMX  propellant  at  Edwards  Air  Force 
Base  (EAFB) ,  and  personnel  there  agreed  to  provide  the  propellant  required 
for  this  program.  After  discussions  with  EAFB  personnel,  we  concluded 
that  the  usual  method  of  constructing  targets  for  Lagrange  gage  experi¬ 
ments  from  machined  pieces  was  not  practical  for  HMX-based  propellants. 
Hence,  an  assembly  based  on  casting  rather  than  machining  the  propellant 
was  designed  to  eliminate  the  complex  machining  and  grooving  operations 
conventionally  used  in  constructing  targets  for  Lagrange  gage  experiments. 
Assemblies  consisting  of  acrylic  tubes  containing  the  stress  gages  and 
ionization  pins  were  constructed  and  sent  to  EAFB  to  be  filled  with 
propellant,  A  gas  gun  experiment  with  Lagrange  stress  gages  embedded  in 
a  propellant  sample  obtained  from  LLNL  was  also  performed  to  test  the 
gage  configurations  in  the  DDT  experiments.  The  pressure  histories 
recorded  in  this  experiment  validated  the  use  of  the  gages  in  the  config¬ 
urations  designed  for  the  DDT  experiments,  DDT  experiments  with  the  filled 

& 

assemblies  obtained  from  EAFB  were  performed  in  steel  and  Lexan  tubes  to: 

•  Record  stress-time  profiles  in  a  propellant  as  it 
underwent  the  transitior  from  deflagration  to 
detonation. 

©  Measure  flame  trajectories  in  propellants  undergoing 
DDT. 

* 

Lexan  is  the  trade  name  for  polycarbonate  manufactured  by  General  Electric 
Co. 
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The  experiments  performed  in  steel  tubes  will  be  presented,  followed  by 
those  performed  in  Lexan  tubes. 

A.  Experiments  in  Steel  Tubes 

The  experimental  arrangement  is  shown  in  Fig.  1.  It  consisted  of  a 
thick-walled  steel  confinement  tube,  end  plugs,  an  igniter  cup,  a  PMMA 
tube  containing  the  propellant  and  stress  gages,  and  bolts  and  nuts  to 
clamp  the  end  plugs  tightly  to  the  confinement  tube.  The  experimental 
arrangement  will  now  be  presented  in  more  detail. 

1.  DDT  Confinement  Tube 

The  confinement  tube  was  made  of  4142  steel.  It  was  28  in. 
long  with  an  inside  diameter  of  3  in.  and  an  outside  diameter  of  7  in. 

The  confinement  tube,  end  plugs,  and  clamping  bolts  were  designed  to 
confine  the  high  gas  pressures,  which  were  generated  by  the  igniter  and 
by  the  burning  propellant,  until  pressures  of  the  order  of  thousands  of 
bars  were  reached.  The  inside  diameter  of  the  tube  was  large  enough  to 
make  any  diameter  effect  on  burning  negligible  and  large  enough  to  exceed 
the  critical  diameter  of  the  propellant.  The  end  plug  at  the  igniter  end 
had  a  hole  through  which  the  igniter  leads  were  passed. 

2.  Igniter 

The  igniter  material  was  a  mixture  of  boron  and  potassium 
nitrate  with  a  laminac-lubersol  binder  obrained  from  Polytechnic 
Specialties  in  Byron,  Georgia,  The  mixture  was  22.7 %  boron,  71.7% 
potassium  nitrate,  and  5.6%  binder.  This  mixture  generates  about 
1600  cal/g  of  material.  About  10  g  of  the  igniter  powder  was  placed 
in  a  Lucite  cup  lined  on  the  bottom  with  a  nichrome  heating  ribbon. 

The  cup  had  an  inside  diameter  of  2.5  in.  and  a  depth  of  0.25  in.  The 
nichrome  ribbon  was  about  11  in.  long,  0,06  in.  wide,  and  0.007  in. 
thick;  it  had  a  resistance  of  about  1  ohm. 
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The  igniter  material  was  thermally  initiated  using  a  28  V,  5  A. ,  dc 
power  supply  to  actuate  the  nichrome  heating  ribbon. 

3.  Propellant  Tube  and  Propellant 

The  Lucite  tube  used  to  contain  the  propellant  inside  the  steel 

DDT  tube  is  shown  in  Fig.  2  as  it  appears  before  the  propellant  is  poured 

into  it.  This  tube  is  a  casting  mold  for  the  propellant  and  a  jig  for 

holding  the  gages  during  casting.  The  main  section  of  the  Lucite  tube  is 

10  in.  long  and  has  an  inside  diameter  of  2  in.  This  section  was  slit 

lengthwise  down  the  center  of  the  cylinder  to  support  the  gages.  Three 

* 

ytterbium  stress  gages  were  glued  to  a  sheet  of  Kapton,  their  leads 
insulated  from  one  another  by  additional  pieces  of  Kapton.  The  resulting 
Kapton  sheet  with  gages  was  then  glued  between  the  two  Lucite  half 
cylinders,  as  shown  in  Fig.  2.  The  stress  gages  were  located  1,  2,  and 
3  in.  from  the  end  of  the  cylinder.  The  gage  leads  were  attached  to  Twinax 
coaxial  cable  and  potted  with  epoxy  inside  a  Lucite  backing  tube  that  had 
an  outside  diameter  of  2.25  in.  (Fig.  2).  The  backing  tube  was  to  delay 
the  arrival  at  the  measuring  stations  of  stress  waves  reflected  from  the 
end  of  the  DDT  tube.  The  entire  assembly,  propellant  and  backing  tube, 
was  22.375  in.  long.  Three  ionization  pins  were  placed  in  the  tube  at 
points  that  were  0.25,  0.5,  and  0.75  in.  from  the  end. 

After  the  stress  gages  and  ionization  pins  were  emplaced,  the 
tubes  were  sent  to  EAFB,  where  the  propellant  was  cast  into  them.  Our 
original  plan  was  no  use  propellant  containing  approximately  50%  HMX. 

However,  because  of  the  propellant  classification  regulations  in  force 
at  the  time,  the  propellant  we  obtained  contained  approximately  12%  HMX. 

A  number  of  1-in. -long  sections  of  Lucite  tube  of  the  same  diameter  as 
the  propellant  tube  were  also  sent  to  EAFB.  These  1-in.  sections,  termed 
add-on  sections,  were  filled  with  either  solid  or  powdered  propellant. 


Kapton  is  a  registered  trademark  of  Du  Pont  for  polyimide  film. 

4» 

Twinak  is  a  two-conductor-shielded  coaxial  cable  manufactured  by 
Trompeter  Electronics. 
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The  powdered  propellant  supplied  for  the  experiments  exhibited  a  rate  of 
pressure  rise  of  2.5  x  10  psi/sec  in  a  quickness  test  performed  by 
personnel  at  EAFB.  Each  of  the  1-in.  add-on  sections  contained  four 
ionization  pins.  The  pins,  0.25-in.  apart,  were  used  to  measure  the 
burning  rate  in  the  propellant.  The  purpose  of  the  add-on  sections  is 
described  in  the  next  section. 


4. 


Experimental  Arrr.i  lgeraent 


Two  1-in, -long  add-on  sections  were  added  in  front  of  the 
10-in.  main  section  of  solid  propellant  in  each  experiment.  We  used 
various  combinations  of  powdered  and  solid  filled  add-on  sections  to 
vary  the  burning  rate  in  the  propellant.  Three  experiments  were 
conducted;  Table  2  gives  the  arrangement  of  the  propellant  add-on  sections. 
The  10-in.  main  section,  in  which  the  stress  gages  were  emplaced,  was 
used  in  all  experiments. 


The  add-on  sections  were  taped  to  the  10- in.  solid-propellant 

section,  and  a  light  pipe  was  epoxied  to  the  outside  of  the  combined 

cylinder,  flush  with  the  igniter  end  of  the  propellant  tube.  The  entire 

propellant  tube — containing  stress  gages,  ionization  pins,  and  light 

pipe— was  then  placed  inside  the  steel  DDT  tube  so  as  to  butt  against 

the  igniter  (see  Fig.  1).  After  alignment,  the  cavity  between  the 

propellant  tube  and  the  inside  diameter  of  the  steel  DDT  cylinder  was 

filled  with  epoxy.  We  used  stycast  low  viscosity  Resin  No.  3020  with 
•k 

catalyst  No.  9.  Because  this  epoxy  has  a  low  exotherm  temperautre 
(V30°C)  it  was  safe  to  use  with  the  propellant.  A  dummy  igniter  was 
used  to  hold  the  components  in  the  proper  positions  during  the  epoxy 
pour.  After  the  epoxy  had  hardened,  the  dummy  igniter  was  removed  and 
the  real  igniter  was  inserted  and  the  end  caps  bolted  together.  The 
bolts,  1.25  in.  in  diameter,  were  torqued  to  1,000  ft-lb  each.  The 
leads  for  the  stress  gages  and  ionization  pins  exited  through  the  end 
cap  opposite  the  igniter  as  shown  in  Fig.  1.  The  space  around  the 
leads  was  filled  with  epoxy. 


From  Emerson  and  Cummings,  Inc.,  Gardena,  CA  90248. 
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Table  2 


ARRANGEMENT  OF  PROPELLANT  SECTIONS 
IN  STEEL  CONFINEMENT 


Experiment  Number 


Propellant  Section 

1 

2 

3 

First  Add-On  Section 

Solida 

,  ,b 

Powdered 

Powdered^ 

Second  Add-On  Section 

Powdered^ 

Powder ed^ 

Solid3 

density  -  1.95  g/cm^. 
^Density  =1.51  g/cm^. 
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5. 


Instrumentation 


The  signals  from  the  stress  gages  and  ionization  pins  were 
recorded  on  oscilloscopes.  The  signal  from  the  light  pipe,  produced 
when  the  burning  in  the  igniter  mixture  reached  the  igniter-propellant 
interface,  was  used  to  trigger  the  oscilloscopes;  the  oscilloscopes 
recorded  the  stress  gage  signals.  The  power  supplies  for  the  stress 
gages  were  turned  on  by  the  "gate  out"  from  the  oscilloscopes  recording 
the  stress  gage  signals.  The  oscilloscopes  that  recorded  the  -ionization 
pins  were  internally  triggered  by  the  signal  coming  from  the  first  pin 
that  closed.  One  oscilloscope,  triggered  by  the  signal  from  the  light 
pipe,  recorded  the  signal  from  the  pins.  This  record,  was  used  to  time- 
correlate  the  stress  gages  and  the  ionization  pins.' 

6.  Results  and  Interpretation  of  Experiments  in  Steel  Tubes 

Table  2  shows  the  arrangement  of  cracked  and  solid  propellant 
in  the  three  experiments  that  were  performed  in  steel  tubes.  The  results 
of  each  experiment  will  be  discussed  and  then  some  general  conclusions 
will  be  given. 

In  the  first  experiment,  in  which  we  had  expected  stresses  in 
the  order  of  40  kbar,  the  stress  gages  recorded  stresses  of  only  about 
5  to  7  kbar.  Such  low  stresses  indicate  that  the  burning  did  not  produce 
high-order  detonation.  Because  the  instrumentation  was  set  for  the 
expected  nigh  stresses,  the  stress  gage  records  obtained  were  not 
optimum  quality.  No  pin  records  were  obtained  from  this  experiment. 

Examination  of  the  recovered  hardware  indicates  that  at  some 
point  during  the  experiment,  venting  occurred  between  the  steel  DDT 
cylinder  and  the  end  cap  at  the  igniter  end.  The  hardware  was  recovered 
still  completely  assembled,  but  the  bolts  had  stretched  under  the  internal 
pressure  generated  by  the  burning  igniter  and  propellant.  After  the 
experiment,  there  was  a  0.125-in.  gap  between  the  igniter  end  cap  and 
the  DDT  cylinder.  Markings  on  the  hardware  indicate  that  the  igniter  - 
cup  was  forced  backward  and  that  the  hot  gases  escaped  between  the 


igniter  cup  and  the  DDT  cylinder.  These  escaping  gases  eroded  the  surfaces 
of  the  cup  and  the  inside  of  the  DDT  cylinder.  Figure  3  shows  the  loose 
fit  between  the  igniter  cup  and  the  steel  cylinder;  the  preshot  tolerance 
of  this  fit  was  0.003  in.  Figures  4(a)  and  (b)  show  that  the  hot  escaping 
gases  also  eroded  the  end  cap  and  the  bolts. 

Although  stronger  bolts  were  used  in  the  second  and  third 
experiments,  venting  still  occurred  and  the  recovered  hardware  looked 
the  same  as  that  in  the  first  experiment.  The  hot  escaping  gases  eroded 
'parts  of  the  end  cap,  the  bolts,  and  the  igniter  cup.  Interestingly, 
we  recovered  pieces  of  Lucite  from  inside  the  confinement  cylinder  that 
did  not  look  as  if  they  had  been  melted. 

Stress  gage  records  were  obtained  from  all  gages  in  the  second  and 
third  experiments.  Figure  5  shows  the  oscilloscope  records  from  Experi¬ 
ment  No.  3  and  Fig.  6  shows  the  corresponding  reduced  stress- time  profiles. 
The  data  in  Fig.  5  and  6  indicate  a  steady  state  of  stress  at  each  gage 
position  about  100  to  200  ysec  after  the  initial  rise.  The  steady  state 
of  stress  decreases  from  gages  1  to  3,  however,  indicating  that  the 
propellant  maintains  a  stress  gradient.  Gage  3  from  experiment  3  recorded 
for  the  longest  time— about  900  ysec. 


The  stress  waves  arrive  at  the  gages  much  sooner  than  the  flame 
front.  In  the  first  experiment,  the  stress  waves  arrived  at  gage  3  (the 
best  record)  about  480  ysec  after  the  flame  reached  the  igniter-propellant 
Interface.  In  the  second  and  third  experiments,  the  stress  waves  arrived  k 

at  the  gages  just  when  the  flame  reached  the  Igniter-propellant  interface. 

Two  reasons  for  this  difference  are  suggested.  In  the  first  experiment,  4 

weaker  bolts  were  used  and  the  venting  may  have  occurred  early,  causing 
a  slower  pressure  buildup.  Also,  in  the  first  experiment  the  first  3 

propellant  add-on  section  after  the  igniter  was  solid  propellant,  whereas 
in  the  second  and  third  experiments  the  first  add-on  propellant  section 
was  powdered  propellant.  Thus,  the  gas  generated  by  the  igniter  might  i 

have  entered  the  propellant  more  readily  in  the  latter  two  experiments, 
resulting  in  the  earlier  generation  of  stress  waves.  1 

'  j 
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(c)  GAGE  3 
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FIGURE  5  OSCILLOSCOPE  RECORDS  FROM  STRESS  GAGES  IN  EXPERIMENT  3 
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FIGURE  6  STRESS  HISTORIES  IN  PROPELLANT  OBTAINED  FROM  GAGE  RECORDS  OF 
EXPERIMENT  3  AND  USED  FOR  LAGRANGE  ANALYSIS 
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Good  quality  ionization  pin  records  were  obtained  for  experi¬ 
ment  No.  3.  Seven  pin  signals  were  recorded;  the  results  are  shown  on 
a  distance-time  plot  in  Fig.  7.  There  is  a  decided  change  in  slope  in 
the  pin  results  at  about  900  ysec.  The  burning  rate  determined  from  the 
pin  data  (ignoring  the  first  pin)  up  to  that  point  was  0.068  mm/ysec. 

After  about  900  ysec,  the  burning  rate  slowed  to  about  0.021  mm/ysec. 

7 .  Lagrange  Analysis 

The  Lagrange  stress  histories  based  on  experimental  data,  shown 
in  Fig.  6,  were  used  in  a  Lagrange  analysis  to  calculate  Lagrange  particle 
velocity  histories,  Lagrange  stress-particle  velocity  paths,  and  Lagrange 
stress-specific  volume  paths  in  compressed  propellant  ahead  of  the  flame 
front.  In  this  Lagrange  analysis,  the  momentum  equation  [equation  (30)] 
is  integrated  with  the  stress  gradient  that  is  estimated  from  the  stress 
records  to  generate  the  particle  velocity  histories  at  the  gage  positions. 
The  continuity  equation  [equation  (29)]  is  then  integrated  with  the 
particle  velocity  gradient  that  is  estimated  from  the  particle  velocity 
histories  to  generate  the  specific  volume  histories  at  the  gage  positions. 
The  elimination  of  time  between  the  stress  histories  and  the  particle 
velocity  histories  gives  the  corresponding  Lagrange  stress-specific  volume 
paths.  The  (u,t)  profiles  calculated  from  the  (a,t)  records  are  shown 
in  Fig.  8,  and  the  corresponding  Lagrange  (a,u)  paths  and  Lagrange  (a,v) 
paths  are  shown  in  Figs.  9  and  10. 

Examination  of  these  profiles  shows  that  the  flame  produces  a 

compression  wave  ahead  of  itself  as  it  propagates  along  the  tube.  The 

stress  at  a  given  gage  position  builds  up  to  a  constant  value  and  a  constant 

Lagrange  pressure  gradient  is  formed  in  the  compression  ahead  of  the  flame. 

A  shock  is  not  formed  in  front  of  the  flame  because,  as  shown  in  Figs.  9 

2  2 

and  10,  the  propellant  satisfies  the  conditions  (3  c/3u  ),<  0  and 

2  2  ” 

(3  o/3v  )^  <  0  in  the  pressure  range  produced  by  the  flame.  The  limited 

buildup  in  pressure  is  attributed  to  the  relatively  small  amount  of  HMX 

in  the  propellant.  Shock  formation  is  expected  at  higher  pressures  in 

this  propellant,  however,  because  these  second  derivatives  become  positive 

as  the  pressure  increases. 
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TIME  FROM  ARRIVAL  OF  FLAME  AT  IGNITER-PROPELLANT  INTERFACE  —  Msec 


LAGRANGIAN  DISTANCE  FROM  IGNITER-PROPELLANT  INTERFACE  —  mm 

MA-6069-10 

FIGURE  7  IONIZATION  PIN  DATA 
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FIGURE  9  LAGRANGE  STRESS-PARTICLE  VELOCITY  PATHS  CALCULATED  FROM  THE 
STRESS  HISTORIES  SHOWN  IN  FIGURE  6 
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FIGURE  10  LAGRANGE  STRESS-VOLUME  PATHS  CALCULATED  FROM  THE  STRESS 
HISTORIES  SHOWN  IN  FIGURE  6 
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Experiments  in  Lexan  Tubes 


The  propellant  KW-6  used  in  these  experiments  contained  50%  HMX. 

The  confinement  was  changed  from  steel  to  Lexan  for  this  propellant  with 
the  higher  percentage  of  HMX  to  make  the  experimental  configuration  much 
lighter  and  more  manageable.  The  experimental  arrangement  is  shown 
schematically  in  Fig.  11.  It  consisted  of  a  thick-walled  confinement 
tube  of  Lexan  polycarbonate,  a  PMMA  tube  containing  ionization  pins  and 
propellant,  an  igniter  cup,  aluminum  end  plugs,  steel  end  plates,  and 
bolts  and  nuts  of  clamp  the  end  caps  and  plates  tightly  to  the  confinement 
tube. 


1.  DDT  Confinement  Tube 

The  confinement  tube  was  made  of  Lexan,  a  polycarbonate.  It 
was  18-in.  long  with  an  outside  diameter  of  6  in.  and  an  inside  diameter 
of  3  in.  The  confinement  tube,  end  plugs,  end  plates,  and  clamping  bolts 
were  designed  to  confine  the  high  gas  pressures,  generated  by  the  igniter 
and  by  the  burning  propellant,  until  pressures  of  the  order  of  thousands 
of  bars  were  reached.  The  inside  diameter  of  the  tube  was  large  enough 
to  make  any  diameter  effect  on  burning  negligible  and  large  enough  to 
exceed  the  critical  diameter  of  the  propellant. 

* 

Two  Viton  o-rings  were  placed  over  the  confinement  tube  at  the 
igniter  end  to  form  a  temporary  gas-tight  seal  between  the  confinement 
tube  and  the  aluminum  end  cap  (see  Fig. 11).  The  high  temperature  per¬ 
formance  of  the  Viton  o-rings  would  seal  against  the  hot  gases  for  the 
life-time  of  the  experiment  (several  seconds) .  Double  o-rings  were  used 
to  extend  the  lifetime  of  the  seal. 

The  end  cap  and  steel  plate  of  the  igniter  end  had  holes  through 
which  the  igniter  leads  and  light  pipe  passed.  At  the  opposite  end,  the 
end  cap  and  steel  plate  had  exit  holes  for  the  ionization  pin  leads  and 
the  stress  gage  leads. 

& 

Viton  -  Trade  name  for  Fluorelastomer  manufactured  by  E.  I.  duPont  de  Nemours 
and  Co .  Inc . 
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FIGURE  11  SCHEMATIC  DRAWING  OF  EXPERIMENTAL  ARRANGEMENT  IN  LEXAN  CONFINEMENT 


A  =  Igniter  Cup  Assembly 

B  =  PMMA  Tube  Containing  Propellant 

C  =  Lexan  DDT  Confinement  Tube 

D  =  End  Caps 

E  =  Light  Pipe 

F  =  High  Strength  Bolts 

G  =  Nuts 

H  =  Steal  End  Plates 
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The  igniter  material  was  a  mixture  of  boron  and  potassium  nitrate 
with  a  laminac-lubersol  binder  obtained  from  Polytechnic  Specialities  in 
Byron,  Georgia.  The  mixture  was- 22.7%  boron,  71.7%  potassium  nitrate, 
and  5.6%  binder.  This  mixture  generates  about  1600  cal/g  of  material. 
Approximately  16  g  of  the  igniter  powder  was  placed  in  a  Lucite  cup  lined 
on  the  bottom  with  nichrome  heating  ribbon.  The  cup  had  an  inside  diameter 
of  2  in.  and  a  depth  of  0.5  in.  The  Lucite  cup  was  placed  in  a  steel  cup 
which  fit  inside  the  Lexan  confinement  tube  as  shown  in  Fig.  11. 

The  nichrome  heater  ribbon  was  about  11-in.  long,  0.06-in.  wide, 
and  0.007-in.  thick,  with  a  resistance  of  about  1  ohm.  The  igniter  leads 
were  brought  out  the  rear  of  the  igniter  assembly  through  the  aluminum 
end  plugs  and  steel  end  plates.  The  igniter  material  was  thermally 
initiated  using  a  28-volt,  5-amp  dc  power  supply  to  heat  the  nichrome 
ribbon. 

3.  Propellant  Tube  and  Propellant 

The  PMfiA  tube  used  to  contain  the  propellant  inside  the  Lexan 
confinement  tube  is  shown  in  Fig,  12.  This  tube  was  to  serve  as  a  casting 
mold  for  the  propellant  and  a  jig  for  holding  the  ytterbium  stress  gages 
during  casting.  The  main  section  of  this  tube  was  2- in.  long  and  the  out¬ 
side  diameter  was  3  in.  This  section  was  slit  lengthwise  down  the  center 
of  the  cylinder  to  support  the  gages. 

* 

Four  ytterbium  stress  gages  were  glued  to  a  sheet  of  Kapton, 
their  leads  were  insulated  from  one  another  by  additional  layers  of  Kapton. 
The  resulting  Kapton  sheet  with  gages  was  then  glued  between  the  two  PMMA 
half  cylinders,  as  shown  in  Fig.  12(a).  The  stress  gages  were  located  1, 

3,  5,  and  7  in,  from  the  end  of  the  cylinder.  The  gage  leads  were  attached 
to  Twinax+  coaxial  cable  and  potted  with  epoxy  inside  the  last  4  in.  of 
the  PMMA  tube.  The  backing  tube  was  to  delay  the  arrival  at  the  gages  of 
stress  waves  reflected  from  the  end  of  the  tube. 

Kapton  is  a  registered  trademark  of  DuPont  for  polyimide  film. 

A. 

Twinax  is  a  two-conductor  shielded  coaxial  cable  manufactured  by  Trompter 

Electronics. 
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The  ytterbium  gages  were  cut  from  0.002-in. -thick  ytterbium 
foil.  Each  gage  was  0.6-in.  long  with  four  legs  (see  Fig.  12(b)).  The 
two  outside  legs  of  each  gage  were  connected  to  the  gage  power  supply. 

The  gage  voltage  was  measured  across  the  two  inside  legs.  The  resistance 
of  each  gage  was  approximately  0.25  £2. 

A  3-in. -long  PMMA  section  was  placed  in  front  of  the  gage  section 
of  the  tube.  This  short  tube  section  containing  ionization  pins  is  shown 
in  place  in  Fig.  12(b).  A  groove  was  cut  in  either  side  of  the  PMMA  pro¬ 
pellant  tube  and  short  tube  in  which  the  ionization  pin  leads  were  placed 
as  shown  in  Fig.  12(b).  By  placing  the  leads  in  grooves,  the  whole 
assembly  could  be  inserted  into  the  confinement  tube. 

Ionization  pins  were  placed  on  either  side  of  the  short  section 
so  that  the  pins  were  1/4-in.  apart  as  shown  in  Figure  13.  The  pins 
were  cemented  in  place  in  steel  collars  with  Eastman  910.  The  lead  ends 
of  the  pins  were  bent  flush  along  the  bottom  of  the  grooves  and  the  leads 
were  epoxied  in  the  grooves  along  the  sides  of  the  short  tube  and  the 
propellant  tube  as  shown  in  Fig.  12(b)  and  Fig.  13. 

After  the  stress  gages  and  ionization  pins  were  emplaced,  the 
tubes  were  sent  to  EAFB,  where  the  propellant  was  cast  into  them.  The 
propellant  KW-6  contained  50%  HMX. 

4.  Experimental  Arrangement 

The  tubes  filled  with  propellant  were  sent  to  us  by  Lt.  Steven 
Cliff  of  EAFB.  The  propellant  was  cast  into  three  of  the  long  PMMA  tubes 
and  two  of  the  short  PMMA  tubes.  The  density  was  2.27  g/cm  .  The  other 
tubes  of  each  kind  were  filled  with  powdered  propellant  packed  to  a  density 
of  65%  of  the  cast  density.  By  using  combinations  of  powdered  and  cast 
propellant,  we  hoped  to  control  the  burning  rate  in  the  propellant  and 
observe  the  effect  on  stress  wave  propagation  and  buildup  at  the  ytterbium 
gage  positions.  Five  experiments  were  conducted  with  the  combinations 
shown  in  Table  3. 


Table  3 


ARRANGEMENT  OF  PROPELLANT  SECTIONS  IN  LEXAN  CONFINEMENT 

Experiment  Number 

1  2  3  4  5 

Cast  Powdered  Powdered  Cast  Powdered 

Cast  Powdered  Powdered  Cast  Cast 
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Propellant  Section 

Short-Pin  Section 
Long-Gage  Section 
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The  short  sections  containing  the  ionization  pins  were  placed 
at  the  appropriate  end  of  the  long  section  containing  the  ytterbium  gages, 
and  the  pin  leads  were  dressed  in  the  grooves  along  the  side  of  the  long 
tubes  as  shown  in  Fig.  12(b).  This  assembly  was  inserted  into  the  Lexan 
confinement  tube.  The  igniter  assembly  was  fitted  in  place,  the  o-rings 
placed  on  the  confinement  tube,  and  the  end  plugs  and  steel  end  plates 
put  in  place  and  bolted  together.  Figure  14(a)  shows  the  components 
arranged  for  assembly  and  Fig.  14(b)  shows  them  assembled  with  2  of  6 
bolts  installed.  The  bolts  were  3/4-in.  diameter  and  were  torqued  to 
500  ft- lb  each. 

5.  Ins  trumen  ta  t ion 

The  signals  from  the  stress  gages  and  ionization  pins  were 
recorded  on  oscilloscopes  and  also  on  a  500-  kilohertz  Fm  tape  recorder. 

The  signal  from  the  light  pipe,  produced  by  the  burning  igniter  mixture, 
triggered  the  electronics.  The  oscilloscopes  recording  the  ionization 
pins  were  triggered  promptly  while  the  oscilloscopes  recording  the  stress 
gage  signals  were  triggered  through  appropriate  delays.  The  power  supplies 
for  the  stress  gages  were  turned  on  by  the  "gate  out"  from  the  oscilloscopes 
recording  the  stress  gage  signals. 

The  tape  recorder  has  essentially  an  "infinite"  time  coverage, 
so  it  was  triggered  promptly  by  the  signal  from  the  light  pipe.  The 
ionization  pins  were  coded  and  recorded  on  two  separate  channels.  Each 
stress  gage  was  recorded  on  a  separate  channel,  The  results  from  the 
tape  recorder  were  used  to  time-correlate  the  stress  gages  and  ionization 
pins. 


6.  Results  and  Interpretation  of  Experiments  in  Lexan  Tubes. 

Table  3  shows  the  arrangement  of  cast  and  powdered  propellant 
in  the  five  experiments  that  were  performed.  Pin  records  were  obtained 
in  all  of  these  experiments  but  ytterbium  gage  records  were  obtained  in 
only  the  first  three.  The  pin  records  were  used  to  calculate  burning 
rates.  The  gage  records  were  used  as  input  in  a  Lagrange  analysis  to 
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integrate  the  equations  of  motion,  equations  (29)  and  (30)  and  determine 
the  states  produced  in  front  of  the  flame.  Values  of  the  particle  velocity 
and  specific  volume  calculated  from  the  gage  records  in  the  Lagrange 
analysis  were,  however,  found  to  be  unrealistic.  The  calculated  particle 
velocities  were  much  too  high  and  the  calculated  specific  volumes  were 
much  too  low.  These  unrealistic  values  were  obtained  because  consecutive 
gage  records  are  separated  by  times  about  a  factor  of  5  greater  than 
those  expected  in  a  one-dimensional  compression  or  rarefaction  wave. 

Because  the  gage  histories  measured  by  the  oscilloscopes  and  the  tape 
recorder  agree,  and  no  such  problem  was  encountered  in  the  DDT  experi¬ 
ments  performed  in  steel  tubes,  we  conclude  that  the  Lexan  tube  did  not 
provide  one-dimensional  confinement  to  the  propellant  in  front  of  the 
flame.  For  this  reason,  we  only  present  here  the  burn  velocities,  the 
maximum  stresses  calculated  from  the  gage  records,  and  the  velocities 
between  gages  in  uncompressed  propellant  calculated  from  the  gage  records. 
These  data  are  given  in  Table  4. 

Examination  of  the  Lexan  tubes  after  the  completion  of  the 
experiments  led  to  the  conclusion  that  while  DDT  did  not  occur  in 
experiments  1,  3,  4,  and  5,  it  may  have  occurred  in  experiment  2  with 
the  highest  calculated  pressure.  Although  the  burning  rate  was  higher 
in  experiment  3  than  in  experiment  2,  venting  occurred  in  experiment  3 
at  the  end  of  the  tube  housing  the  gage-leads. 

We  conclude  that  Lexan  does  not  provide  the  confinement  required 
for  a  successful  Lagrange  gage  study  of  DDT,  and  we  recommend  the  use  of 
steel-confining  tubes  in  any  such  study  that  may  be  undertaken  in  the 
future. 
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Table  4 


EXPERIMENTAL  RESULTS  FOR  LEXAN  CONFINEMENT 


Burn 

Maximum 

Wave  . 

Velocity 

Pressure 

Velocity 

(mm/ys) 

(kbar) 

(mm/ys) 

0.066 

Experiment  1 

0.045 

1.9 

0.5 

Experiment  2 

0.12 

>35 

0.22 

Experiment  3 

0.18 

4+ 

0.24 

Experiment  4 

0.07 

a 

- 

Experiment  5 

0.06 

_a 

- 

From  arrival  lines. 

+Gage  Stretched. 

£ 

No  gage  records  were  obtained. 
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Model  solutions  are  presented  for  the  deflagration  to  shock  transition  (DST) 
process  associated  with  the  propagation  of  an  accelerating  flame  in  propel¬ 
lant  contained  in  a  closed  tube.  The  flame  is  treated  as  a  reactive  discon¬ 
tinuity,  the  flow  in  front  of  the  flame  as  a  simple  compression  wave,  and 
the  density  gradient  in  the  flow  behind  the  flame  is  assumed  to  be  zero. 

The  solutions  are  based  on  expressions  for  the  burning  rate  and  the  accel¬ 
eration  of  the  flame  and  are  valid  until  the  time  the  reaction  becomes  com¬ 
plete  in  the  flame.  Consideration  of  the  evolution  of  the  shock  preceding 
the  flame  after  this  time  led  to  a  postulate  for  assessing  the  probable 
occurrence  of  the  deflagration  to  detonation  transition  (DDT).  Examination 
of  particular  solutions  in  the  light  of  this  postulate  showed  that  DDT  de¬ 
pends  strongly  on  the  burning  rate  and  will  probably  not  occur  unless  the 
burning  rate  parameter  exceeds  a  critical  value. 


a  shock  wave  and  the  properties  of  this  shock 
satisfy  the  conditions  required  for  the  occur¬ 
rence  of  SDT.  The  condition  for  shock  formation 
is  thus  not  a  criticality  condition  for  DDT,  and 
the  shock  must  be  characterized  to  assess  the 
probability  of  the  occurrence  of  DDT,  It  is  for 
this  reason  that  our  present  phenomenological 
treatment  of  DDT  is  restricted  primarily  to  an 
analytic  treatment  of  DST. 


TREATMENT  OF  DDT 


INTRODUCTION 

The  deflagration  to  detonation  transition 
(DDT)  in  a  propellant  or  explosive  is  said  to 
occur  when  the  propagation  of  a  flame  results  in 
the  formation  of  a  detonation  wave  in  front  of 
the  flame.  Significant  events  envisaged  in  the 
evolution  of  DDT  are  as  follows: 

•  The  formation  of  an  accelerating  flame 
by  an  unsteady  burning  process. 

•  The  formation  of  a  shock  wave  in  unburnt 
material  by  the  compressive  action  of  the 
accelerating  flame  in  the  deflagration  to 
shock  transition  (DST) . 

•  The  formation  of  a  detonation  wave  in 
material  ahead  of  the  flame  by  the  shock 
to  detonation  transition  (SDT) . 

Vie  are  concerned  in  the  present  paper  with 
the  occurrence  of  DDT  in  a  burning  propellant  con¬ 
tained  in  a  closed  tube.  Our  approach  to  this 
problem  is  based  on  the  premise  that  it  is  valid 
to  model  DDT  by  the  above  chain  of  events  when 
the  burning  process  is  subject  to  the  confinement 
provided  by  a  closed  tube.  In  this  case,  DDT  will 
not  occur  unless  the  accelerating  flame  produces 


This  work  was  supported  by  the  Air  Force  Office 
of  Scientific  Research  under  Contract  F49620-77- 
C-0039. 

65 


Our  approach  is  based  on  the  premise  that 
the  dynamics  of  the  flow  behind  and  in  front  of 
an  accelerating  flame  must  be  considered  to  char¬ 
acterize  the  shock  formed  in  front  of  the  flame 
and  assess  the  probability  of  DDT.  In  this  case, 
an  analytical  treatment  of  DDT  must  model  the 
accelerating  flame,  and  the  unsteady  flows  on  both 
sides  of  it.  The  model  for  these  unsteady  pro¬ 
cesses  used  in  the  present  paper  is  based  on  the 
model  formulated  by  Jones  [1]  to  investigate  the 
DDT  process  in  gases.  We  treat  the  flame  as  a 
reactive  discontinuity,  the  flow  in  front  of  the 
flame  as  a  simple  compression  wave,  and  we  assume 
that  the  density  gradient  is  zero  in  the  flow 
behind  the  flame.  The  one-dimensional  differen¬ 
tial  equations  governing  this  type  of  flow  in  a 
closed  tube  can  be  solved  when  a  relationship  for 
the  burning  rate  is  known  and  the  acceleration  of 
the  flame  is  prescribed.  The  validity  of  such 
solutions  must  be  questioned  ss  the  flow  develops 
however  because  the  equation  used  to  prescribe  the 


; 
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acceleration  of  the  flame  may  become  incompatible 
with  the  Hugoniot  equation  expressing  the  balance 
of  energy  across  the  flame.  The  solution  pre¬ 
sented  here  is  only  valid  until  the  time  (tj>) 
the  reaction  becomes  complete  in  flame  because 
the  Hugoniot  equation  shows  that  the  extent  of 
reaction  increases  in  the  flame  as  the  flow  de¬ 
velops.  Particular  solutions  were  constructed 
and  examined  over  the  time  intervals  defined  by 
the  reaction  times  to  determine  the  times  taken 
for  shock  formation  in  front  of  the  flame,  and  the 
position  and  strength  of  these  shocks  at  later 
times.  The  strength  of  the  shock  when  all  the 
chemical  energy  is  liberated  in  the  flame  is  then 
used  to  relate  the  probability  of  DDT  to  the 
burning  rate.  Our  treatment  is  phenomenological 
but  supercedes  previous  thermo-hydrodynamic  treat¬ 
ments  [2-4]  by  making  the  states  attained  across 
the  flame  compatible  with  the  differential  equa¬ 
tions  governing  the  conservation  of  mass  and  mo¬ 
mentum  in  the  flow  behind  the  flame  as  well  as  with 
the  rear-boundary  particle  velocity  condition  in 
the  closed  tube. 


SOLUTION  FOR  DST 

We  let  the  superscript  f  and  F  denote  res¬ 
pectively  the  flame  discontinuity  and  its  propa¬ 
gation  velocity,  and  we  let  the  subscripts  +  and 
-  denote  quantities  in  front  of  and  behind  the 
flame.  States  across  the  flame  are  related  by 
the  Rankine-Hugoniot  jump  conditions  written  as 


u  +  c  ■  c  +  2u 
+  +  o  + 


We  denote  the  extent  and  heat  of  reaction  by  X 
and  q,  assume  that  the  reaction  products  are 
polytropic  with  an  index  k  -  3,  and  write  the 
equation  of  state  for  material  behind  the  flame 
as 


e  -  Aq  +  • 
o  ^ 


Ap  p  v  (p  ) 

.  too  _  (1  _  x) 


2  0 
,c  p  -2p  N 

+  (1  -  X)(-2-e - 2>(V+  (p_)  -  Vq) 


where  v+(p_)  is  written  to  denote  that  the  volume 
of  the  propellant  is  related  to  the  pressure  by 
the  same  equation  on  both  sides  of  the  flame  and 
vP  denotes  the  specific  volume  of  the  detonation 
products  in  their  standard  state.  The  combination 
of  Eqs.  (3),  (6),  and  (8)  allows  us  to  write  the 
extent  of  reaction  in  the  flame  as 

Xf  -  AEf/Qf  (9) 

where  Qf  -  [2q  +  pQ(v£  -  Vq)  +  (c^pQ  -  3pQ)  x 

(v+(pf)  -  vQ)]  9a 

and  AEf  -  [(p_  +  p+)(v_  -  v+)  +  (pv)_  -  (pv)+  + 

(copo  "  3V  (v+(P>v0>]f  % 


P*(F  -  UJ)  =  pf (F  -  uf)  (1) 

P+  -  pf  P+(F  -  "  ul)  (2) 

2(e*  -  ef+)  =  (p|  +  p^)(v*  -  v*)  (3) 

where  P,  u,  p,  and  e  denote  density,  particle 
velocity,  pressure  and  specific  energy,  and  v  •• 
I/O  denotes  the  specific  volume.  We  let  c  dt 
the  sound  speed  and  the  subscript  o  denote  the 
initial  state  in  the  propellant,  and  we  assume 
that  the  propellant  is  governed  by  a  Murnaghan 
equation  with  an  index  n  »  3.  The  equations  for 
p+,  p+,  e  ,  and  the  slope  of  a  forward  facing  C+ 
characteristics  ii+  +  c+)in  the  simple  isentropic 
compression  wave  can  then  be  written  as 


2  [  /  \  3 

cp  /  u  \ 

0  0  /,  +1 

+  ”  [\  " 


1  +  po  (4) 


e.  ■  e  + 
+  o 


<P+  -  P0)v+ 
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with  f  written  as  a  superscript  on  the  bracket 
for  notational  convenience. 

We  let  x  and  t  denote  Eulerian  distance  and 
time  so  that  F  -  dxf/dt.  In  contrast  to  Jones  [1], 
we  assume  that  the  burning  rate  is  related  to  the 
particle  velocity  and  the  pressure  at  the  flame 
front  by  the  equation 

F  «  u^  +  A(p^  -  PQ)  (10) 

where  A  is  a  parameter.  The  assumption  that  the 
pressure  exponent  in  Eq.  (10)  is  1  is  not  neces¬ 
sary  but  is  made  here  for  convenience. 

We  prescribe  the  acceleration  of  the  flame 
by  the  assumption  that  the  particle  velocity  at 
the  flame  front  is  governed  by  a  power  law  of  the 
form 


with  Oq  and  m  considered  as  parameters. 

We  are  now  in  a  position  to  construct  the 
solution  for  the  flow  produced  by  our  accelera¬ 
ting  flame.  We  write  the  differential  equations 
for  the  flow  in  front  of  the  flame  aa 


(p  c**  -  2p  ) 

,00  0/  s  ,  y  V 

+  - - -  Cv+  -  v0)  (6) 


3p+  3p+  3u,, 

IT  +  u+  IT  +  P+  IT 
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+  u 


+  3x 


+  P 


ri4\  *t 

+  ydpy+  3x 


(13) 


to  satisfy  the  isentropic  condition,  and  we  write 
the  differential  equation  for  the  flow  behind  the 
flame  as 


3p_ 


3u 


3t  +  p-  3x 


(14) 


3u_ 

3 r 


+  u 


3u 

3x 


+  P 


3p_ 

3x 


(15) 


to  satisfy  the  assumption  that  3p_/3x  «  0  behind 
the  flame.  It  is  convenient  to  construct  the 
solution  in  terms  of  u£  rather  than  t.  The  pro¬ 
cedure  is  as  follows:  (1)  calculate  the.  flame 
trajectory,  (2)  calculate  the  time  for  shock  for¬ 
mation  t8  in  front  of  the  flame,  and  (3)  calcu¬ 
late  the  density  behind  the  flame  to  determine 
the  burnt  state,  the  u_(x,t)  and  p_(x,t)  fields 
behind  the  flame,  and  the  reaction  time  t^ . 


We  will  now  carry  out  these  steps  in  detail: 
(1)  Calculation  of  the  Flame  Trajectory 


We  first  combine  Eqs.  (4)  and  (10)  to 
derive  the  following  equation  for  the  flame  speed 
in  terms  of  uj, 

'  2  .  ,.3-1 


+  3 


(16) 


where  B  =  A  cqPq/3.  The  exponent  of  the  pressure 
in  Eq.  (10)  was  assumed  to  be  1  in  order  to  ob¬ 
tain  this  simple  expression  for  F.  We  combine 
Eq.  (16)  with  the  identity 


dt 


(17) 


du. 


to  obtain  a  differential  equation  for  x  and  in¬ 
tegrate  it  with 


1-m 


(18) 


and  mC  «  (c^)  ^B(c0)^1  derived  from  Eq.  (11) 
to  obtain  the  following  expression  for  x^, 


1+m 


[ (3  +  co/B) 


1  + 


(19) 
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The  equation-  for  the  position  of  the  flame  as  a 
function  of  time  is  simply  obtained  by  combining 
Eqs.  (19)  and  (11). 


(2)  Shock  Formation  in  Front  of  the  Flame 


A  shock  is  formed  in  the  compression 
wave  ahead  of  an  accelerating  flame  when  the  pro¬ 
pellant  satisfies  the  condition  (d2p/dp2)+  >  0  be¬ 
cause  the  C+  characteristics  emanating  from  the 
flame  front  form  an  envelope.  Since  the  discon¬ 
tinuity  occurs  in  the  flow  where  3x/3u+  <=  32x/3u2 
“  0,  application  of  these  conditions  to  the  par¬ 
ticle  velocity  field  determines  the  position  x8 
and  the  time  ts  of  its  formation.  When  Eqs.  (12) 
and  (13)  are  integrated  by  the  method  of  charac¬ 
teristics  from  the  flame  front,  the  particle 
velocity  field  in  the  compression  wave  can  be 
written  as 


(c  +  2uf)t  +  f(u^) 
o  +  + 


(20) 


and  x  and 


tf  de- 


with  f(uF)  -  xf  -  (c0  +  2u*)tf  „„  „  . . 

fined  in  terms  of  u£  by  Eqs.  (19)  and  (11).  For¬ 
mulating  the  conditions  3x/3u^  «  32x/3(uf)2  =  0 
with  Eqs.  (20),  (19),  and  (11)  leads  after  some 
manipulation  to  the  equations 


tf- 


^)-(co  +  2uJ)Vtf 


*du, 


(21) 


/2(m  +  l)uf  f  muf 
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dF 
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f  -  Co/x 
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8(u+) 


where 


d2tf 


A!  M 

d(u+) 


-^-^andtf 
mu+  du+ 


f  dt 


+  ,  f 
du+ 


2  f  f  2 

When  m  »  1,  and  d  t  /d(u+)  »  0  everywhere  in  the 

flow,  the  equation  for  the  envelope  of  the  char¬ 


acteristics  written  es 


cot  +  xf  +  (2u^  -  F(u£))tf 


(23) 


shows  that  the  shock  is  formed  on  the  first  char¬ 
acteristic  where  ■  uj  •  0  *  Setting 

t^  «  h£  ■*  u£  “  F(u£)  ■  0  in  Eq.  (21)  then  gives 
the  time  the  shock  is  formed  on  the  first  char¬ 
acteristic  as 
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When  i  =f  1,  the  condition 


2(m  +  l)u 


(m  -  1) 


l+  f  mu+ 
±  +  ?(V-7^r 


dF 


duf 


(25) 


mist  be  solved  to  find  the  particle  velocity 
(u|.)8,  and  the  characteristic  associated  with  the 
discontinuity.  The  time  the  shock  is  formed  ts 
can  then  be  calculated  from  Eq.  (21)  and  the  po¬ 
sition  where  it  is  formed  can  be  calculated  from 
Eq.  (20).  Attention  in  this  paper  is  focused  on 
the  case  when  m  »  1. 

The  evolution  of  the  shock  after  t?  is 
determined  by  applying  a  graphical  method  [5]  to 
the  region  of  the  simple  wave  where  the  particle 
velocity  becomes  double  valued.  In  this  graphi¬ 
cal  method  based  on  a  weak  shock  approximation, 
the  position  of  the  discontinuity  in  tne  (u+  -  x) 
plane  is  determined  by  the  condition  that  the 
ureas  between  the  (u+  -  x)  profile  and  the  dis¬ 
continuity  are  the  same  on  both  sides  of  the  dis¬ 
continuity. 


(3)  Calculation  of  States  Behind  the  Flame 

The  calculation  of  statec  behind  the  flame 
is  based  on  the  assumption  that  (3p_/3x)  “  0. 

In  this  case,  the  density  p_  can  be  readily  cal¬ 
culated  from  an  equation  expressing  the  conserva¬ 
tion  of  mass.  Since  p_  *  p_,  the  calculated  den¬ 
sity  can  be  used  to  calculate  u£  and  p£  from  the 
jump  conditions,  Eqs.  (1)  and  (2)  ,  and  these 
values  can  be  used  to  complete  the  characteriza¬ 
tion  of  the  u  (x,t)  and  p_  (x,t)  fieldsbehind  the 
flame.  The  following  expression  for  p_, 

f 


f 


f  f.  dt  .  f 

P+  <*-«+>  7T  du+ 

du+ 


(26) 


is  obtained  as  an  integral  expression  for  the 
conservation  of  mass  by  equating  the  mass  behind 
the  flame  with  the  mass  that  passed  through  the 
flame  front.  The  mass  behind  the  flame  can  be 
written  simply  a r.  o_x^  because  p_  can  be  taken 

r*f 

outside  the  integral  I  p_dx  when  (3p_/3x)  «  0. 
o 

When  the  integration  in  Eq.  (26)  is  performed, 

P  can  be  written  in  terms  of  u£  as 


B  C  c  m 

P_(u+)  «  p, 


(14m)  (l+2m)\c 


(143m) 


f  f 

with  x  (u+)  given  by  Eq 


I  . 

2 

f\ 

3 

I  — 

+  1 . 

± 

\  C 

\  0 

(144.., 

lCoi 

(27) 


(19) .  The  combination 
of  Eqs.  (27)  and  (19)  gives  the  equation  for  the 
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in  terras  of  uj.  When  p_/p 
is  known,  the  following  equations  obtained  from 
Eqs. 


(28) 


density  ratio  p_/p 
'  iown,  the  folio 
(1),  (2)  and  (10), 


f  f  f  /  f  \2  /P4/P0  \ 

p.  - p+  - e+  (A p+)  ter-y 


(29) 


can  be  used  to  calculate  i/  and  as  a  function 
of  u^. 

The  assumption  that  p_  »  p  (t)  is  also  used 
to  construct  the  equations  for  the  particle  velo¬ 
city  and  pressure  behind  the  flame.  It  follows 
from  Eq.  (14)  that  the  particle  velocity  gradient 
behind  the  flame  is  only  a  function  of  time  and 
can  be  written  as 


'  3u 
,  9x 


n(t) 


(30) 


Integrating  Eq.  (30)  and  making  use  of  the  boun¬ 
dary  condition  u  *  0  at  x  »  0  gives  the  equation 
for  the  particle  velocity  field  behind  the  flame 
as 


ri(t)x 


(31) 


where  h(t)  -  u_(t)/x  (t)  and  0  <  x  <  x^ .  Inte¬ 
gration  of  Hq.  (15)  with  the  expression  for  the 
acceleration  derived  from  Eq.  (31)  then  gives 
the  equation  for  the  pressure  field  behind  the 
flame  as 

p_»  pf(t)  +  Y  |n2  +  £^((xf(t)2-  X2)  (32) 

When  p1  and  v1  have  been  calculated  as  functions 
of  u£,\he  time  tR  for  reaction  to  go  to  comple¬ 
tion  in  the  flame  can  be  calculated  from  Eq.  (11) 
with  the  value  of  u£  that  gives  X*  «  1  in  Eq.  (9). 


PARTICULAR  SOLUTIONS  FOR  DST  AND  THEIR  APPLICA¬ 
TION  TO  DDT 


:o  the  following  set  or  parameters  tp  c  ,  q,  vq, 
n,  am,  A) .  The  solutions  presented  here  nave  the 

2  mm/we,  q  »  54.69  kbar 


Particular  solutions  for  DST  are  obtained 
from  the  previous  equations  by  assigning  values 
to  the  following  set  of  parameters  (po,  co,  q, 
m, 

values  R0  «  1.6  g/cc,  c0  .  .  .  , 

cc/g,  vP  •  8.4  x  102cc/g,  m  •  1,  a  »  oij .  *  2  x  10' 2 
mm/us2,  and  values  of  the  burning  parameter  A  in 
the  range  9.0  X  10“3-  1.2  x  10"2  mm/us  kbar.  The 
value  q  is  based  on  an  assumed  Chapman- Jouquet 
pressure  in  the  propellant  of  350  kbar*  Equations 
for  the  flow  variables  were  evaluated  for  specified 
values  of  u£.  Equations  (U)  and  (19)  were  used  to 
calculate  the  flame  path,  Eqs.  (4)  through  (7)  and 
Eq.  (20)  were  used  to  calculate  states  in  the 
compression  wave,  Eqs.  (19),  (27),  (28),  (29), 


] 

I 


$  . 

}\  . 


and  Eq.  (9)  were  used  to  calculate  values  of  the 
flow  variables  in  the  burnt  state  and  the  density 
behind  the  flame,  and  Eq.  (31)  was  used  to  calcu¬ 
late  values  of  the  particle  velocity  behind  the 
flame.  The  value  of  was  increased  in  these 
calculations  until  the  reaction  coordinate  satis¬ 
fied  the  condition  »  1  to  determine  the  reaction 
tR  and  define  the  time  limit  imposed  on  the  solu¬ 
tion  by  the  energy  equation. 

Plots  of  the  particle  velocity-distance  prr 
files  calculated  with  the  following  values  of  A, 
1.2  x  10'2  mm/us  kbar,  1.1  x  10-2  mm/ys  kbar,  and 
9.0  x  10-3  mm/ys  kbar  are  shown  in  Figures  1,  2, 
and  3  to  demonstrate  the  dependence  of  DST  on 
this  burning  rate  parameter.  The  flame  path  and 
the  states  connected  by  the  flame  with  A  =  1.2  x 
1CT2  mm/ys  kbar  are  given  in  Table  1.  For  conve¬ 
nience  in  presentation.  Figure  1  is  shown  on  this 
page,  but  Figures  2  and  3  and  the  table  are  shown 
on  the  next  page. 


c 

shown  at  20  ys,  35  ys,  and  52  ys.  Because  \  »  1 

in  the  discontinuity  shown  as  F^F",  the  reaction 
time  t^  «  52  ys  for  this  flow.  The  particle  velo-  ' 
city  profiles  FFl  and  FF"  behind  the  flame  are 
shown  as  straight  lines  to  satisfy  Eq.  (30).  Be¬ 
cause  m  -  1,  the  shock  is  formed  on  the  wave  front 
at  S  where  the  particle  velocity  gradient  becomes 
infinite  after  the  time  t8  «  50  ps  given  by  Eq. 

(24) .  The  flow  therefore  satisfies  the  condition 
t8  <  tp  because  the  shock  is  formed  before  the 
reaction  goes  to  completion  in  the  flame  discon¬ 
tinuity.  Our  consideration  of  the  shock  after  t^, 
when  the  prescription  of  the  flame  acceleration 
is  no  longer  valid,  is  based  on  the  assumption  that 
the  flame  can  be  treated  as  a  constant  velocity  or 
decelerating  piston.  In  either  case,  the  informa¬ 
tion  that  the  flame  has  undergone  a  change  in  its 
motion  will  be  propagated  along  the  C+  character¬ 
istic  emanating  from  the  flame  front  when  t  «=  tR. 
Points  in  the  (u+  -  x)  profile  F"C  in  Figure  1 
will  thus  continue  to  propagate  along  C+  character¬ 


Figure  1.  Flame  path,  shock  path,  and  particle  velocity-distance  profiles  for 

the  flow  produced  by  an  accelerating  flame  with  A  »  1.2  x  10'2  mm/us  kbar. 


The  particle  velocity-distance  profiles  shown 
in  Figure  1  will  now  be  considered  in  more  detail 
to  gain  an  insight  into  DST  and  formulate  a  pos¬ 
tulate  for  assessing  the  probability  of  DDT. 

Figure  1  shows  the  formation  and  development  of 
the  shock  produced  by  the  accelerating  flame  when 
A  »  1.2  x  10“  mm/ys  kbar.  The  flame  path  is  shown 
as  FFjF^  and  the  shock  path  as  SS'.  The  posi¬ 
tions  of  the  flame  discontinuity  and  the  particle 
velocity  profiles  on  both  aides  of  the  flame  are 


istics  with  velocities  given  by  Ea.  (7).  The  shock 
particle  velocity  and  the  position  of  the  shock 
after  a  time  t  >tg  »  50  ys  are  defined  by  equalizing 
the  areas  in  the  double  valued  regions  of  the 
(u+  -  x)  profile  as  discussed  earlier.  The  (u+  -  x) 
profile  S'C'  at  t  ■  60  ys  in  Figure  1  for  example 
gives  0.2  mm/ys  for  the  shock  particle  velocity  and 
121  mm  for  the  position  of  the  shock.  The  corre¬ 
sponding  value  of  the  shock  pressure  at  t  «  60  ys 
is  7.1  kbar.  Application  of  this  graphical  method 
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Figure  2: 


Parcicle  velocity-distance 
profiles  for  the  flow  when 
the  shock  is  formed  and  the 
reaction  becomes  complete  in 
the  accelerating  flame  with 
A  *»  1.1  x  i0“2  mm/us  kbar. 


Figure  3:  Particle  velocity-distance 
profile  for  the  flow  when 
the  reaction  is  complete  in 
the  accelerating  flame  with 
A  =  9.0  x  10“3  mm/ps  kbar. 
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Table  1.  Calculated  values  of  the  position  of  the  flame 
and  the  states  attained  across  the  accelerating 
flame  discontinuity  with  A  *  l.  2  x  10"^  smt/us  kbar. 


to  define  the  shock  as  the  flow  develops  after  tR 
leads  to  the  conclusion  that  the  final  shock  pressure 
will  be  the  pressure  p+  »  53.6  kbar  attained  at  the 
flame  front  when  C  -  tR.  Consequently,  in  applying 
these  solutions  for  an  accelerating  flame  to  DDT, 
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we  postulate  that  the  pressure  attained  at  the  flame 
front  when  the  reaction  proceeds  to  completion  in 
the  flame  can  be  used  to  assess  the  probability  of 
the  occurrence  of  DDT. 


.  V  h-.v  ■■  -  ,  - 

•  ..f  V 


Particle  velocity-distance  profiles  for  flows  •  A  burning  rate  governing  by  the  square 

when  A  =  1.1  x  10"2mm/us  kbar  and  A  ■=  9.0  x  10-3  of  the  pressure  is  nit  necessary  for 

mm/ vis  kbar  are  shown  respectively  in  Figures  2  the  onset  of  DDT  in  a  closed  tube, 

and  3.  The  flame  trajectory  and  the  particle 

velocity  profile  at  tR  is  shown  in  both  figures.  The  model  solutions  constructed  in  this  paper 

and  it  is  clear  from  these  profiles  that  tR  <  ts  for  the  flow  associated  with  an  accelerating  re- 

for  both  these  values  of  A.  Figure  2  also  shows  active  flame  discontinuity  provide  guidelines  for 

the  profile  in  the  compression  wave  at  50  ys  when  assessing  the  probability  of  the  occurrence  of  DDT 
the  shock  is  formed.  According  to  our  postulate,  in  a  burning  propellant  contained  in  a  closed  tube, 

the  shock  pressures  for  assessing  the  probability  Questions  regarding  the  applicability  of  these 

of  DDT  are  respectively  '37  kbar  and  -25  kbar  3  solutions  to  real  systems  must  be  considered  how- 

when  A  =  1.1  x  10“  mtn/ps  kbar  and  A  =  9.0  x  10~  ever  because  the  extent  of  reaction  in  the  flame 

mm/us  kbar.  discontinuity  increases  as  the  flow  develops.  It 

is  my  conjecture  that  the  assumptions  made  to 

The  shock  pressures  for  assessing  the  proba-  treat  the  flow  may  be  valid  in  a  damaged  or  cast 

biUty  of  DDT  is  our  three  flows  are  therefore  burning  composite  propellant  when  the  components 

54  kbar,  -37  kbar,  and  -25  kbar.  While  DDT  is  have  markedly  different  reaction  rates, 

expected  to  occur  when  the  shock  pressure  is  54 

i  kbar  or  37  kbar,  it  may  not  occur  when  the  shock  A  more  fundamental  theoretical  study  of  the 

pressure  is  25  kbar.  The  dependence  of  the  pres-  flow  produced  by  an  unsteady  flame  must  be  under- 
sure  at  the  flame  front  on  A  when  t  -  tg  shows  taken  to  obtain  a  basic  understanding  of  DDT. 

that  the  probability  of  the  occurrence  of  DDT  Questions  regarding  the  treatment  of  a  flame  as 

depends  strongly  on  the  burning  rate,  and  leads  a  reactive  discontinuity  and  the  application  of 

to  the  conclusion  that  DDT  will  probably  not  occur  a  steady-state  burning  rate  expression  to  an  un¬ 
in  a  closed  tube  when  A  <  0.9  x  10"  ^  mm/ps  kbar.  steady  flame  must  be  addressed  in  this  study,  and 

a  more  complete  treatment  of  the  flow  must  be  in¬ 
corporated  into  a  model  for  predicting  the  accel- 
RESULTS  AND  CONCLUSIONS  eration  of  the  flame. 


Solutions  describing  the  one-dimensional  un¬ 
steady  flow  produced  by  an  accelerating  flame 
propagating  in  a  propellant  or  explosive  contained 
in  a  closed  tube  were  constructed  to  model  DST  and 
gain  an  insight  into  DDT.  The  flame  was  treated 
as  a  reactive  discontinuity,  the  flow  in  front  of 
the  flame  was  treated  as  a  simple  compression 
wave,  and  the  density  gradient  was  assumed  to  be 
zero  in  the  flow  behind  the  flame.  The  descrip¬ 
tion  of  the  flow  is  phenomenological  because  the 
burning  rate,  the  acceleration  of  the  flame,  and 
the  density  gradient  behind  the  flame  must  be 
prescribed  to  construct  the  solution,  but  it 
supercedes  thermodynamic  treatments  by  accounting 
for  the  flow  behind  the  flame  as  well  as  the  flow 
in  front  of  it.  A  graphical  method  was  applied 
to  the  particle  velocity  profiles  in  the  compres¬ 
sion  wave  to  determine  the  development  of  the 
shock  ahead  of  the  flame.  Particular  solutions 
were  constructed  over  the  time  interval  0  <  t  < 
tg  defined  by  the  time  tg  taken  for  the  reaction 
to  proceed  to  completion  in  the  flame  disconti¬ 
nuity.  It  is  not  worthwhile  to  consider  detailed 
properties  of  these  solutions  after  the  reaction 
time  tg  because  the  equation  used  to  prescribe 
the  acceleration  of  the  flame  is  not  valid  after 
this  time.  The  postulate  that  the  pressure  at¬ 
tained  at  the  flame  front  when  t  ■  tg  can  be  used 
to  assess  the  probability  of  DDT  was  based  on  the 
assumption  that  the  flame  can  be  treated  as  a 
constant  velocity  or  decelerating  piston  after 
tg,  Examination  of  particular  solutions  in  the 
light  of  this  postulate  led  to  the  following 
conclusions: 

•  DDT  depends  strongly  on  the  burning 
rate  and  will  probably  not  occur 
unless  the  burning  rate  parameter 
exceeds  a  critical  value. 
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